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This  Semiannual  Report  also  presents  statistics  from  operation  of  the  Intelligent  Monitor¬ 
ing  System  (IMS).  The  IMS  has  been  operated  in  an  experimental  mode  using  NORESS 
and  ARCESS  data,  and  the  performance  has  been  very  satisfactory.  Preparations  have 
been  made  for  installing  a  new  version  of  the  IMS  that  accepts  data  from  an  arbitrary  num¬ 
ber  of  arrays  and  single  3-component  stations. 

The  NORSAR  Detection  Processing  system  has  been  operated  throughout  the  period  with 
an  average  uptime  of  99.3%  as  compared  to  98.6%  for  the  previous  reporting  period.  A 
total  of  21 38  seismic  events  have  been  reported  in  the  NORSAR  monthly  seismic  bulletin. 
The  petformance  of  the  continuous  alarm  system  and  the  automatic  bulletin  transfer  by 
telex  to  AFTAC  has  been  satisfactory.  A  system  for  direct  retrieval  of  NORSAR  wave¬ 
form  data  through  an  X.25  connection  has  been  implemented,  and  has  been  tested  success¬ 
fully  for  acquiring  such  data  by  AFTAC.  Processing  of  requests  for  full  NORSAR/ 
NORESS  data  on  magnetic  tapes  has  progressed  according  to  established  schedules. 

On-line  detection  processing  and  data  recording  at  the  NORSAR  Data  Processing  Center 
(NDPC)  of  NORESS,  ARCESS,  FINES  A  and  GERESS  data  have  been  conducted 
throughout  the  period.  Data  from  the  two  stations  in  Poland  have  been  recorded  and  pro¬ 
cessed  in  an  expeiimental  mode.  Monthly  processing  statistics  for  the  arrays  as  well  as 
results  of  the  IMS  analysis  for  the  reporting  period  are  given. 

There  have  been  n  lodifications  made  to  the  NORSAR  data  acquisition  system.  The  pro¬ 
cess  of  evaluating  and  testing  technical  options  for  refurbishment  of  the  array  is  continu¬ 
ing. 

Maintenance  activities  in  the  period  comprise  preventive/conective  maintenance  in  con¬ 
nection  with  all  the  NORSAR  subanays,  NORESS  and  ARCESS.  In  addition,  the  mainte¬ 
nance  center  has  been  inv()lved  with  occasional  maintenance  of  equipment  for  FINESA 
and  work  in  connection  with  the  two  stations  in  Poland.  Other  activities  have  involved 
testing  of  the  NORSAR  communicati()ns  systems,  and  field  studies  at  sites  in  Spitsbet^en 
and  the  Kola  Peninsula. 

Summaries  of  seven  scientific  contributions  are  presented  in  Chapter  7  of  this  report. 

In  Section  7.1.  the  application  of  statistical  parameter  estimation  theory  to  .seismic  holog¬ 
raphy  algorithms  is  described.  The  theoretical  framework  for  this  application  is  devel¬ 
oped,  and  the  advantages  of  the  statistical  algorithms  are  analyzed  by  computer  modelling 
experiments.  Data  from  the  six  sensors  of  NORSAR  subairay  01 A  for  a  Hindu  Kush 
earthquake  (22  February  1972  m^,  =  5.6)  have  been  used  in  this  experiment  together  with 
simulated  scattered  waveforms.  The  statistical  algorithm  is  found  to  be  superior  to  the 
conventional  seismic  holography  algorithm  in  this  case,  as  it  provides  better  resolution 
and  suppresses  ''leakage”  effects  in  the  imaging  process. 
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Section  7.2  describes  a  new  maximum  likelihood  scheme  for  3-component  slowness  esti¬ 
mation.  This  new  technique  incorporates  the  response  of  the  crust  beneath  the  receiver  in 
the  signal  model.  The  paper  presents  examples  of  applications  both  for  a  single  3C  station 


and  an  array  of  four  3C  stations,  and  demonstrates  the  utility  of  the  method  for  slowness 
analysis  of  P-wave  recordings.  Plans  are  being  made  to  extend  the  technique  to  analysis  of 
SV  and  SH  waves. 

Section  7.3  follows  up  our  previous  studies  on  using  adaptive  statistically  optimal  algo¬ 
rithms  for  on-line  event  detection.  It  is  found  that  regional  P-phase  detection  at  typical  fre¬ 
quencies  (2.5-4.5  Hz)  requires  noise  adaptation  for  8  azimuths  and  2  velocities  in  order  to 
achieve  less  than  3  dB  worst  case  missteering  loss.  Preliminary  results  of  applying  a  one¬ 
dimensional  statistically  optimal  detector  (ESTDET)  are  presented  and  compared  with  the 
traditional  STA/LTA  algorithm.  The  ESTDET  method  is  found  to  be  more  sensitive,  at  a 
similar  level  of  false  alarm  rates. 

In  section  7.4  we  describe  the  current  status  and  planned  extensions  of  the  regional  array 
network  in  northern  Europe.  A  new  digital  3-component  station  in  Apatity,  USSR,  is 
described,  and  preliminary  analysis  shows  that  it  is  very  valuable  for  near-field  studies  of 
mining  explosions  in  the  Khibiny  Massif.  Noise  surveys  and  plans  for  a  small-aperture 
array  at  the  Spitsbergen  archipelago  are  described.  This  array  is  expected  to  become  oper¬ 
ational  in  1992. 

Section  7.5  presents  an  application  of  statistically  optimal  adaptive  processing  of  array 
data  in  the  context  of  separating  interfering  events.  Besides  the  obvious  importance  of  this 
application  in  detecting  possible  explosions  in  the  coda  of  large  earthquakes,  it  also  serves 
to  improve  the  capability  of  the  "Threshold  Monitoring”  technique.  This  preliminary 
study  indicates  that  with  the  appropriate  modifications,  the  adaptive  processing  may  sig¬ 
nificantly  improve  the  detection  threshold  during  event  coda  conditions. 


In  Section  7.6.  a  case  study  is  presented  on  regional  detection  performance  during  the 
recently  conducted  GSETT-2  experiment.  Using  the  bulletin  of  the  Seismological  Insti¬ 
tute,  University  of  Helsinki,  as  a  reference,  it  is  shown  that  the  90%  detection  threshold  for 
Fennoscandia/NW  Russia  is  close  to  2.5  in  teims  of  the  duration  magnitudes  used  in  that 
bulletin.  This  shows  that  in  a  region  with  dense  coverage  of  high  quality  arrays  such  as 
Fennoscandia,  an  excellent  monitoring  capability  may  be  achieved.  However,  the  exten¬ 
sion  of  these  results  to  different  geological  environments  in  other  parts  of  the  world  is 
unproven. 


Section  7.7  describes  the  contemporary  seismicity  of  the  NW  part  of  the  USSR.  The  paper 

is  based  on  recordings  from  the  regional  station  network  operated  by  the  Kola  Science 

Centre  since  1982.  Selected  seismic  zones  are  described  in  some  detail;  the  Kola- 

Finnmark  zone,  the  Khibiny  Massif,  NE  Archangelsk  and  a  seismic  zone  in  Novaya 

Zemlya.  The  Khibiny  Massif  zone  is  of  particular  interest,  since  it  appears  that  a  great 

increase  of  induced  seismicity  has  taken  place  in  recent  years,  following  the  increased 
•  .1.  For 

mining  activity  in  the  area.  ;  . . 
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1  Summary 

This  Semiannual  Technical  Summary  describes  the  operation,  maintenance  and  research 
activities  at  the  Norwegian  Seismic  Array  (NORSAR),  the  Norwegian  Regional  Seismic 
Array  (NORESS)  and  the  Arctic  Region^  Seismic  Array  (ARCESS)  for  the  period  1 
April  ~  30  September  1991.  Statistics  are  also  presented  for  additional  seismic  stations, 
which  through  cooperative  agreements  with  institutions  in  the  host  countries  provide  con¬ 
tinuous  data  to  the  NORSAR  Data  Processing  Center  (NPDC).  These  stations  comprise 
the  Finnish  Experimental  Seismic  Array  (FINESA),  the  German  Experimental  Seismic 
Array  (GERESS),  and  two  3-component  stations  in  Poland:  Ksiaz  and  Stary  Folwark. 

This  Semiannual  Report  also  presents  statistics  from  operation  of  the  Intelligent  Monitor¬ 
ing  System  (IMS).  The  IMS  has  been  operated  in  an  experimental  mode  using  NORESS 
and  ARCESS  data,  and  the  performance  has  been  very  satisfactory.  F*reparations  have 
been  made  for  installing  a  new  version  of  the  IMS  that  accepts  data  from  an  arbitrary  num¬ 
ber  of  arrays  and  single  3-component  stations. 

The  NORSAR  Detection  Processing  system  has  been  operated  throughout  the  period  with 
an  average  uptime  of  99.3%  as  compared  to  98.6%  for  the  previous  reporting  period.  A 
total  of  21 38  seismic  events  have  been  reported  in  the  NORSAR  monthly  seismic  bulletin. 
The  performance  of  the  continuous  alarm  system  and  the  automatic  bulletin  transfer  by 
telex  to  AFTAC  has  been  satisfactory.  A  system  for  direct  retrieval  of  NORSAR  wave¬ 
form  data  through  an  X.25  connection  has  been  implemented,  and  has  been  tested  success¬ 
fully  for  acquiring  such  data  by  AFTAC.  Processing  of  requests  for  full  NORSAR/ 
NORESS  data  on  magnetic  tapes  has  progressed  according  to  established  schedules. 

On-line  detection  processing  and  data  recording  at  the  NORSAR  Data  Processing  Center 
(NDPC)  of  NORESS,  ARCESS,  FINESA  and  GERESS  data  have  been  conducted 
throughout  the  period.  Data  from  the  two  stations  in  Poland  have  been  recorded  and  pro¬ 
cessed  in  an  experimental  mode.  Monthly  processing  statistics  for  the  arrays  as  well  as 
results  of  the  IMS  analysis  for  the  reporting  period  are  given. 

There  have  been  no  modifications  made  to  the  NORSAR  data  acquisition  system.  The  pro¬ 
cess  of  evaluating  and  testing  technical  options  for  refurbishment  of  the  array  is  continu¬ 
ing. 

Maintenance  activities  in  the  period  comprise  preventive/corrective  maintenance  in  con¬ 
nection  with  all  the  NORSAR  suban-ays,  NORESS  and  ARCESS.  In  addition,  the  mainte¬ 
nance  center  has  been  involved  with  occasional  maintenance  of  equipment  for  FINESA 
and  work  in  connection  with  the  two  stations  in  Poland.  Other  activities  have  involved 
testing  of  the  NORSAR  communications  systems,  and  field  studies  at  sites  in  Spitsbeigen 
and  the  Kola  Penin.sula. 

Summaries  of  seven  scientific  contributions  are  presented  in  Chapter  7  of  this  report. 

In  Section  7. 1 ,  the  application  of  statistical  parameter  estimation  theory  to  seismic  holog¬ 
raphy  algorithms  is  described.  The  theoretical  framework  for  this  application  is  devel- 
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oped,  and  the  advantages  of  the  statistical  algorithms  are  analyzed  by  computer  modelling 
experiments.  Data  from  the  six  sensors  of  NORSAR  subarray  01 A  for  a  Hindu  Kush 
ea^quake  (22  February  1972,  m^  =  5.6)  have  been  used  in  this  experiment  together  with 
simulated  scattered  waveforms.  The  statistical  algorithm  is  found  to  be  superior  to  the 
conventional  seismic  holography  algorithm  in  this  case,  as  it  provides  better  resolution 
and  suppresses  “leakage”  effects  in  the  imaging  process. 

Section  7.2  describes  a  new  maximum  likelihood  scheme  for  3-component  slowness  esti¬ 
mation.  This  new  technique  incorporates  the  response  of  the  crust  beneath  the  receiver  in 
the  signal  model.  The  paper  presents  examples  of  applications  both  for  a  single  3C  station 
and  an  array  of  four  3C  stations,  and  demonstrates  the  utility  of  the  method  for  slowness 
analysis  of  P-wave  recordings.  Plans  are  being  made  to  extend  the  technique  to  analysis  of 
SV  and  SH  waves. 

Section  7.3  follows  up  our  previous  studies  on  using  adaptive  statistically  optimal  algo¬ 
rithms  for  on-line  event  detection.  It  is  found  that  regional  P-phase  detection  at  typical  fre¬ 
quencies  (2.5-4.5  Hz)  requires  noise  adaptation  for  8  azimuths  and  2  velocities  in  order  to 
achieve  less  than  3  dB  worst  case  missteering  loss.  Preliminary  results  of  applying  a  one¬ 
dimensional  statistically  optimal  detector  (ESTDET)  are  presented  and  compared  with  the 
traditional  STA/LTA  algorithm.  The  ESTDET  method  is  found  to  be  more  sensitive,  at  a 
similar  level  of  false  alarm  rates. 

In  section  7.4  we  describe  the  current  status  and  planned  extensions  of  the  regional  array 
network  in  northern  Europe.  A  new  digital  3-component  station  in  Apatity,  USSR,  is 
described,  and  preliminary  analysis  shows  that  it  is  very  valuable  for  near-field  studies  of 
mining  explosions  in  the  Khibiny  Massif.  Noise  surveys  and  plans  for  a  small -aperture 
array  at  the  Spitsbei^gen  archipelago  are  described.  This  array  is  expected  to  become  oper¬ 
ational  in  1992. 

Section  7.5  presents  an  application  of  statistically  optimal  adaptive  processing  of  array 
data  in  the  context  of  separating  interfering  events.  Besides  the  obvious  importance  of  this 
application  in  detecting  possible  explosions  in  the  coda  of  large  earthquakes,  it  also  serves 
to  improve  the  capability  of  the  “Threshold  Monitoring”  technique.  This  preliminary 
study  indicates  that  with  the  appropriate  modifications,  the  adaptive  processing  may  sig¬ 
nificantly  improve  the  detection  threshold  during  event  coda  conditions. 

In  Section  7.6,  a  case  study  is  presented  on  regional  detection  performance  during  the 
recently  conducted  GSETT-2  experiment  Using  the  bulletin  of  the  Seismological  Insti¬ 
tute,  University  of  Helsinki,  as  a  reference,  it  is  shown  that  the  90%  detection  threshold  for 
Fennoscandia/NW  Russia  is  close  to  2.5  in  terms  of  the  duration  magnitudes  used  in  that 
bulletin.  This  shows  that  in  a  region  with  dense  coverage  of  high  quality  arrays  such  as 
Fennoscandia,  an  excellent  monitoring  capability  may  be  achieved.  However,  the  exten¬ 
sion  of  these  results  to  different  geological  environments  in  other  parts  of  the  world  is 
unproven. 

Section  7.7  describes  the  contemporary  seismicity  of  the  NW  part  of  the  USSR.  The  paper 
is  based  on  recordings  from  the  regional  station  network  operated  by  the  Kola  Science 
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Centre  since  1982.  Selected  seismic  zones  are  described  in  some  detail;  the  Kola- 
Finnmark  zone,  the  Khibiny  Massif,  NE  Archangelsk  and  a  seismic  zone  in  Novaya 
Zemlya.  The  Khibiny  Massif  zone  is  of  particular  interest,  since  it  appears  that  a  great 
increase  of  induced  seismicity  has  taken  place  in  recent  years,  following  the  increased 
mining  activity  in  the  area. 


NORSAR  Sci.  Rq>.  1-91/92 


November  1991 


2  NORSAR  Operation 

2.1  Detection  Processor  (DP)  operation 

There  have  been  37  breaks  in  the  otherwise  continuous  operation  of  the  NORSAR  online 
system  within  the  current  6-month  reporting  interval.  The  uptime  percentage  for  he 
period  is  99.3%  as  compared  to  98.6%  for  the  previous  period. 

Fig.  2.1.1  and  the  accompanying  Table  2.1.1  both  show  the  daily  DP  downtime  for  the 
days  between  1  April  and  30  September  1991.  The  monthly  recording  times  and  percent¬ 
ages  are  given  in  Table  2.1.2. 


The  breaks  can  be  grouped  as  follows: 

a)  Hardware  failure  8 

b)  Stops  related  to  program  work  or  error  0 

c)  Hardware  maintenance  stops  0 

d)  Power  jumps  and  breaks  0 

e)  TOD  error  correction  0 

0  Communication  lines  29 


The  total  downtime  for  the  period  was  28  hours  and  56  minutes.  The  mean-tirae-between- 
failures  (MTBF)  was  5.0  days,  as  compared  to  4.1  for  the  previous  period. 


J.  Torstveit 


Fig.  2.1.1.  Detection  Processor  uptime  for  April  (top).  May  (middle)  and  June  (bottom) 
1991. 
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fig.  2.1.1.  Detection  Processor  uj^me  for  July  (top),  August  (middle)  and  Septembn* 
(bottom)  1991. 
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Date 

Time 

Cause 

18  May 

1220  -  1832 

Hardware  failure 

18  May 

1900  -  1932 

Hardware  failure 

18  May 

2038  -  2145 

Hardware  failure 

13  Aug 

1043  -  1248 

Hardware  failure 

24  Aug 

0704  -  0807 

Hardware  failure 

04  Sep 

0821  -  0848 

Hardware  failure 

06  Sep 

0952  -  1400 

Hardware  failure 

10  Sep 

1834  - 

Line  failure 

11  Sep 

-  0453 

16  Sep 

0848  -  1021 

Line  failure 

24  Sep 

1916  -  1959 

Hardware  failure 

Table  2.1.1.  The  major  downtimes  in 

the  period  1  April  -  30  September  1991. 

Month  DP  Uptime 

DP  Uptime 

No.  of 

No.  of  Days 

DP  MTBF* 

Hours 

% 

DP  Breaks 

with  Breaks 

(days) 

Apr  91 

719.83 

99.98 

5 

3 

5.0 

May  91 

735.85 

98.90 

6 

4 

4.4 

Jun  91 

719.73 

99.967 

7 

7 

3.7 

Jul91 

743.50 

99.93 

11 

7 

2.6 

Aug  91 

740.87 

99.58 

2 

2 

10.3 

Sep  91 

703.28 

97.66 

6 

6 

4.2 

99.34 

37 

29 

5.0 

’Mean-timc-bctwecn-failurcs  =  total  uplinK/no.  of  up  intervals. 


Table  2.1.2.  Online  system  performance,  1  April  -  30  September  1991. 


NORSAR  Sci.  Rep.  1-91/92 


November  1991 


2.2  Array  communications 

General 

Table  2.2.1  reflects  the  performance  of  the  communications  system  throughout  the  report¬ 
ing  period.  The  most  common  events  which  have  affected  the  NORSAR  array  have  been: 
damaged  communications  cable,  NTA  transmission  system  failure,  synchronization  prob¬ 
lems,  bad  connections,  attenuator  fault,  lightning  and  low  line  level. 


Detailed  Summary 

April  (weeks  14-17),  1-28.4.91 

03C,  06C  were  affected  (week  15),  03C  by  synchronization  problems,  and  06C  probably 
by  a  short  line  outage.  The  other  systems  performed  satisfactorily. 

May  (weeks  18-22),  29.4-2.6.91 

Most  systems  were  affected  in  May.  It  started  with  06C  (week  19)  which  lost  synchroniza¬ 
tion.  but  after  a  Modcomp  restart  the  system  resumed  operation.  All  systems  were  affected 
18  May  (week  20).  Observations  made  indicated  irregularities  related  to  the  NTA  trans¬ 
mission  systems. 

Attempts  to  start  the  lines  failed.  01  A,  03C  and  04C  were  then  turned  off  (i.e.,  no  com¬ 
mand  words  transmitted  to  the  subarrays).  A  Modcomp  restart  started  the  lines  and  resyn¬ 
chronized  the  remaining  systems:  01 B,  02B,  02C  and  06C. 

21  May  (week  21)  the  three  subarrays  (OlA,  03C  and  04C)  resumed  operation  after  new 
attempts.  04C  failed  again  the  same  day  (21  May),  and  we  were  not  able  to  restart  the  sys¬ 
tem. 

23  May  (week  21 )  the  subarray  was  visited.  Loop  tests  were  run  in  B-  and  C-Ioop  (local 
and  digital).  C-loop  was  accepted,  but  B-loop  failed,  which  was  reasonable  since  part  of 
the  04C  data  path  between  Hamar  and  Berger  Central  now  is  on  PCM  equipment. 
Although  the  communications  systems  were  inoperative  before  the  test  started,  the  C-loop 
functioned  well  when  activated.  An  answer  to  this  may  be  that  activities  in  the  CTV  prior 
to  the  test  reset  the  system. 

31  May  (week  22).  Data  from  OlA  were  not  accepted.  The  subarray  was  masked. 

June  (weeks  23-26),  3-30.6.91 

We  had  synchronization  problems  regarding  02B,  02C  and  03C  (weeks  24,  25  and  26). 
Apparently,  02B  and  later  03C  cau.sed  the  trouble.  We  believe  02C  lost  synchronization  in 
connection  with  the  Modcomp 's  heavy  engagement  in  trying  to  synchronize  02B. 
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13  June  02B  was  turned  off  here  at  NORSAR  (i.e.,  isolated  from  the  Modcomp  I/F  and  the 
remainder  of  the  02B  communications  path). 

27  June  02B  was  visited  and  tested  in  the  B/C-loop,  with  the  test  equipment  hooked  up  to 
the  NDPC  modem.  The  test  did  not  reveal  reduced  system  quality,  and  was  therefore 
switched  on  again. 

02C  was  turned  off  14  June,  and  after  a  Modcomp  restart,  the  system  resumed  operation. 
Tlie  03C  communications  system  was  not  quite  reliable,  as  it  was  difficult  to  resynchro¬ 
nize  after  loss  of  synch.  Also  after  a  small  TOD  correction  the  system  lost  synch. 

July  (weeks  27-31),  1.7-4.8.91 

Only  two  systems  were  affected  this  month:  01 A  and  03C.  The  latter  was  affected  week 
27  due  to  a  defective  attenuator  located  in  Elverum.  Week  29  we  lost  contact  with  OlA, 
caused  by  a  defective  communications  cable  in  the  OlA  area.  03C  again  lost  synch,  week 
30.  A  Modcomp  restart  resynchronized  the  system. 

August  (weeks  32-35),  5.8-1.9.91 

Three  systems  were  affected  this  period:  02B  week  35  by  a  bad  connection  located  in  the 
NTA  premises  Lillestrom;  02C  by  a  broken  communications  cable  weeks  32, 33  and  34  in 
the  Mesnalia  area;  and  04C  week  33,  probably  in  connection  with  an  NTA  transmission 
equipment  outage. 

September  (weeks  36-39),  2-29.9.91 

In  September  all  systems  were  affected,  either  individually  or  all  simultaneously.  02B  out¬ 
age  continued  through  week  35,  and  the  subarray  resumed  operation  2  September.  OlA 
was  down  between  6  and  9  September  because  of  NTA  cable  work  in  the  area.  01 B  was 
down  a  few  hours  on  6  September. 

In  connection  with  service  on  a  4381  disk  6  September,  the  system  was  taken  down  for 
approximately  1  hour.  In  that  period  no  data  acquisition  took  place. 

10  September  (after  3  p.m.)  lines  1  and  3  (OlA,  02B)  stopped.  During  the  period  the  Mod¬ 
comp  tried  to  resynch.  OlA,  02B;  the  remaining  subarrays  lost  synch  also.  02B  was  proba¬ 
bly  the  main  cau.se  of  the  outages,  due  to  line  degradation  and  loss  of  power  in  connection 
with  lightning  in  the  area. 

1 1  September  at  5  a.m.  the  Modcomp  was  restarted  with  02B  turned  off  (i.e.,  no  transmis¬ 
sions  of  ICTWs  toward  the  subarray).  An  almost  identical  situation  occurred  a  few  hours 
later  (after  8  a.m.).  01 A  and  01 B  stopped,  followed  by  03C.  The  remaining  subarrays 
“went  down”  approximately  1 1  a.m.  (02B  was  already  down).  Two  Modcomp  restarts 
were  necessary  in  order  to  restore  operation.  01 B  communications  cable  was  seriously 
damaged  by  a  grounding  spear,  and  the  communications  path  toward  the  subarray  was 
therefore  broken. 
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12  September  approx.  0445  a.m.  03C  resynchronized.  1 3  September  01 B  communications 
cable  was  repaired. 

Finally,  on  16  September  0845  a.m.  we  again  lost  data  from  all  subarrays  (02B  was  still 
down).  According  to  NTA/Hamar  the  reason  was  an  inoperative  transmission  path 
between  Hamar  and  Lillestr0m.  Three  Modcomp  restarts  had  to  be  carried  out  before  all 
lines  started  (-02B),  approximately  10  a.m.  19  September  02B  restored  after  SLEM  reset. 

O.A.  Hansen 
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2.3  Event  Detection  operation 

In  Table  2.3.1  some  monthly  statistics  of  the  Detection  and  Event  Processor  operation  are 
given.  The  table  lists  the  total  number  of  detections  (DPX)  triggered  by  the  on-line  detec¬ 
tor,  the  total  number  of  detections  processed  by  the  automatic  event  processor  (EPX)  and 
the  total  number  of  events  accepted  after  analyst  review  (teleseismic  phases,  core  phases 
and  total). 


Total 

DPX 

Total 

EPX 

P-phases 

Accepted  Phases 

Core  Phases  Sum 

Daily 

Apr  91 

10050 

1316 

337 

51 

388 

12.9 

May  91 

8400 

1050 

270 

48 

318 

10.3 

Jun91 

7225 

1213 

412 

51 

463 

15.4 

Jul91 

8225 

1336 

313 

70 

383 

12.4 

Aug  91 

9350 

1275 

259 

82 

341 

11.0 

Sep  91 

9775 

1123 

200 

45 

245 

8.2 

1791 

347 

2138 

11.7 

Table  2J.1.  Detection  and  Event  Proces.sor  statistics.  1  April  -  30  September  1991 . 


NORSAR  Detections 

The  number  of  detections  (phases)  reported  by  the  NORSAR  detector  during  day  091 
1991,  through  day  273  1991,  was  52,231,  giving  an  average  of  289  detections  per  pro- 
ces.sed  day  (181  days  processed).  The  distribution  of  these  detections  is  given  in  Table 
2.3.2. 

B.  Paulsen  Gammelby 
J.  Fyen 
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3  Operation  of  regional  arrays 

3.1  Recording  of  NORESS  data  at  NDPC,  Kjeller 


Table  3.1.1  lists  the  main  outage  times  and  reasons,  and  as  can  be  seen  the  main  reason  for 
the  outage  is  damage  to  the  Hub  caused  by  a  thunderstorm. 

The  average  recording  time  was  89.89%  as  compared  to  99.2%  during  the  previous  report- 


ing  period. 

Date 

Time 

Cause 

07  Apr 

2106 

-  2138 

Transmission  line  failure 

13  May 

2142 

-  2158 

Transmission  line  failure 

05  Jun 

0505 

-  1055 

Software  work 

05  Jun 

1105 

-  1127 

Software  work 

18  Jun 

1735- 

HUB  failure  due  to  thunderstorm 

25  Jun 

-  1158 

15  Jul 

1122 

-  1254 

Transmission  line  failure 

15  Jul 

1320 

-  1346 

Transmission  line  failure 

26  Jul 

2048 

-  2052 

Transmission  line  failure 

27  Jul 

0550 

-  0605 

Transmission  line  failure 

31  Jul 

2230 

- 

HUB  failure 

02  Aug 

-  0720 

16  Aug 

1445 

- 

HUB  failure  due  to  thunderstorm 

26  Aug 

-  1400 

28  Aug 

0555 

-  0559 

Transmission  line  maintenance 

28  Aug 

0622 

-  0626 

Transmission  line  maintenance 

28  Aug 

0759 

-  0825 

Transmission  line  maintenance 

28  Aug 

0837 

-  0856 

Transmission  line  maintenance 

10  Sep 

1228 

-  1245 

Transmission  line  failure 

17  Sep 

0624 

-  0632 

Transmission  line  failure 

18  Sep 

1118 

-  1126 

Transmission  line  maintenance 

Table  3.1.1.  Interruptions  in  NORESS  recordings  at  NDPC,  1  April  -  30  September  1991. 
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Monthly  uptimes  for  the  NORESS  on-line  data  recording  task,  taking  into  account  all  fac¬ 
tors  (field  installations,  transmission  line,  data  center  operation)  affecting  this  task  were  as 
follows: 


April 

99.91 

May 

99.% 

June 

76.58 

July 

99.50 

August 

63.50 

September 

99.91 

Fig.  3.1.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
NORESS  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period. 


November  1991 


Fig.  3.1.1.  (cont.)  NORESS  data  recording  uptime  for  July  (top),  August  (middle)  and 
September  (bottom)  1991. 
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3.2  Recording  of  ARCESS  data  at  NDPC,  Kjeller 

The  main  reasons  causing  most  of  the  ARCESS  outage  in  the  period  are:  service  on  the 
satellite  link,  and  hardware  failure  at  NDPC. 

The  average  recording  time  was  97.32%  as  compared  to  98.4%  for  the  previous  reporting 
period.. 


Date  Time  Cause 


03  Apr 

0933 

-  0946 

Service  on  the  antenna 

04  Apr 

1511 

-  1544 

Hardware  failure  NDPC 

04  Apr 

1612 

-  1645 

Hardware  failure  NDPC 

04  Apr 

1715 

-  1745 

Hardware  failure  NDPC 

04  Apr 

1810 

-  1840 

Hardware  failure  NDPC 

04  Apr 

1901 

-  1911 

Hardware  failure  NDPC 

04  Apr 

1924 

- 

Hardware  failure  NDPC 

05  Apr 

-  0555 

05  Apr 

0620 

-  0629 

Hardware  failure  NDPC 

05  Apr 

0643 

-  1019 

Hardware  failure  NDPC 

10  Apr 

0807 

-  0820 

Service  satellite  link 

10  Apr 

0826 

-  0838 

Service  satellite  link 

10  Apr 

0856 

-  1009 

Service  satellite  link 

10  Apr 

1043 

-  1140 

Service  .satellite  link 

10  Apr 

1359 

-  1407 

Service  satellite  link 

1 1  Apr 

0702 

-  0717 

Service  satellite  link 

17  Apr 

0132 

-  1535 

Power  break  HUB 

22  Apr 

2204 

- 

Software  failure 

23  Apr 

-  1013 

01  Jun 

0516 

-  0534 

Hardware  failure 

06  Jun 

0735 

-  0819 

Hardware  failure 

17  Jun 

0953 

-  1626 

Service  HUB 

17  Jun 

1658 

-  1709 

Service  HUB 

17  Jun 

1725 

-  1759 

Service  HUB 

17  Jun 

2040 

-  2102 

Service  HUB 

18  Jun 

0839 

-  0904 

Service  HUB 

18  Jun 

1156 

-  1313 

Service  HUB 

19  Jun 

0445 

-  0453 

Software  service 

19  Jun 

0459 

-  1200 

Software  service 

18  Jul 

1035 

-  1043 

Unknown 

25  Jul 

2037 

-  2108 

Unknown 

26  Jul 

0629 

-  0757 

Hardware  service 

26  Jul 

0804 

-  1226 

Hardware  service 

12  Aug 

1807 

-  1830 

Software  failure 
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14  Aug 

0625 

-  0824 

14  Aug 

0900 

-  1134 

14  Aug 

1146 

-  1202 

15  Aug 

0653 

-  0702 

15  Aug 

0858 

-  0908 

16  Aug 

0157 

- 

17  Aug 

-  1217 

24  Aug 

1916 

-  1924 

24  Aug 

1931 

-  1838 

15  Sep 

1321 

-  2004 

Service  satellite  link  at  HUB 
Service  satellite  link  at  HUB 
Service  satellite  link  at  HUB 
Service  satellite  link  at  HUB 
Service  satellite  link  at  HUB 
Satellite  link  failure 

Satellite  link  failure 
Satellite  link  failure 
Power  break  HUB 


Table  3.2.1.  The  main  interruptions  in  ARCESS  recordings  at  NDPC,  1  April  -  30  Sep¬ 
tember  1991. 


Monthly  uptimes  for  the  ARCESS  on-line  data  recording  task,  taking  into  account  all  fac¬ 
tors  (field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were 
as  follows: 


April 

93.62% 

May 

99.99% 

June 

97.56% 

July 

99.12% 

August 

94.59% 

September 

99.06% 

Fig.  3.2.1 .  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
ARCESS  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period. 
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Fig.  3.2.1.  ARCESS  data  recording  uptime  for  July  (top),  August  (middle)  and  September 
(bottom)  1991. 
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33  Recording  of  FINESA  data  at  NDPC,  Kjeller 

The  average  recording  time  was  66.2%  as  compared  to  98.9  for  the  previous  period.  As 
can  be  seen  from  Table  3.3.1  below,  the  main  reason  for  the  downtime  is  transmission  line 
failure  and  HUB  failure. 


Date  Time  Cause 


05  Apr 

1030 

-  1129 

Transmission  line  failure 

16  Apr 

1149 

-  1252 

Service  at  HUB 

27  Apr 

1241 

-  1407 

Hardware  failure  NDPC 

28  Apr 

2015 

-  2113 

Hardware  failure  NDF*C 

03  May 

0053 

-  0556 

Hardware  failure  NDPC 

05  May 

1515 

-  1552 

Hardware  failure  NDPC 

06  May 

0259 

-  0356 

Hardware  failure  NDPC 

06  May 

1031 

-  1052 

Hardware  failure  NDPC 

07  May 

1407 

-  1502 

Hardware  failure  NDPC 

08  May 

0756 

-  0812 

Hardware  service  NDPC 

08  May 

1046 

-  1302 

Hardware  service  NDPC 

14  May 

2047 

-  2127 

Transmission  line  failure 

!5  May 

1509 

- 

Transmission  line  failure  due  to  thunderstorm 

17  May 

-  0755 

Transmission  line  failure  due  to  thunderstorm 

17  May 

1355 

- 

Transmission  line  failure  due  to  thunderstorm 

20  May 

-  1008 

Transmission  line  failure  due  to  thunderstorm 

22  May 

1249 

- 

Transmission  line  failure 

23  May 

-  1027 

Transmission  line  failure 

25  May 

1140 

- 

Bad  data  due  to  numerous  gaps 

28  May 

-  1110 

Bad  data  due  to  numerous  gaps 

30  May 

0841 

-  1201 

Software  upgrade 

03  Jun 

0642 

-  1308 

Software  upgrade 

08  Jun 

0319 

-  0658 

Transmission  line  failure 

10  Jun 

0941 

- 

Transmission  line  service 

12  Jun 

-  0900 

Transmission  line  service 

12  Jun 

0900 

- 

Transmission  line  failure 

26  Jun 

-  0822 

Transmission  line  failure 

29  Jun 

0001 

- 

Transmis.sion  line  failure 

17  Jul 

- 

0620 

Transmission  line  failure 

17  Jul 

1618 

- 

Transmission  line  failure 

18  Jul 

-  0934 

Transmission  line  failure 

18  Jul 

1033 

- 

Transmission  line  failure 

22  Jul 

-  0950 

Transmission  line  failure 

29  Jul 

1100 

-  1408 

Transmission  line  failure 

03  Aug 

1625 

- 

Transmission  line  failure 

05  Aug 

-  1221 

Transmission  line  failure 
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05  Aug 
06  Aug 

1316 

-  0630 

09  Aug 

0452 

-  0937 

09  Aug 

12  Aug 

2243 

-  1148 

12  Aug 

16  Aug 

1439 

-  0706 

16  Aug 
22  Aug 

1048 

-  0955 

27  Aug 

1908 

-  1937 

12  Sep 

0912 

-  1055 

13  Sep 

0627 

-  0853 

Transmission  line  failure 
Transmission  line  failure 
Transmission  line  failure 
Transmission  line  failure 
Transmission  line  failure 
Transmission  line  failure 
Transmission  line  failure 
Transmission  line  failure 
Transmission  line  failure 
Transmission  line  failure 
Transmission  line  test 
Service  HUB 


Table  3  J.l.  The  main  interruptions  in  FINESA  recordings  at  NDPC,  1  April  -  30  Septem¬ 
ber  1991. 


Monthly  uptimes  for  the  FINESA  on-line  data  recording  task,  taking  into  account  all  fac¬ 
tors  (field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were 
as  follows: 


April 

99.29% 

May 

75.64% 

June 

38.76% 

July 

32.09% 

August 

51.86% 

September 

99.39% 

Fig.  3.4. 1  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
FINESA  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period. 


J.  Torstveit 


Fig.  3J.1.  (cont)  FINESA  data  recording  uptime  for  July  (top),  August  (middle)  and  Sep¬ 
tember  (bottom)  1991. 
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3.4  Event  detection  operation 

This  section  reports  results  from  simple  one-array  automatic  processing  using  signal  pro¬ 
cessing  recipes  and  ‘ronap’  recipes  for  the  ep  program  (NORSAR  Sci.  Rep.  No  2-88,89). 

IMS  results  are  reported  in  section  3.5. 


NORESS  detections 

The  number  of  detections  (phases)  reported  during  day  091  1991,  through  day  273  1991, 
was  31,966,  giving  an  average  of  193  detections  per  processed  day  (166  days  processed). 

Table  3.4.1  shows  daily  and  hourly  distribution  of  detections  for  NORESS. 

Events  automatically  located  by  NORESS 

During  days  091  1991,  through  273  1991, 2218  local  and  regional  events  were  located  by 
NORESS,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  1 3.4  events  per  processed  day  (166  days  processed).  63  %  of  these  events  are  within 
300  km,  and  88  %  of  these  events  are  within  1000  km. 


ARCESS  detections 

The  number  of  detections  (phases)  reported  during  day  091  1991,  through  day  273  1991, 
was  84,791,  giving  an  average  of  463  detections  per  processed  day  (183  days  processed). 

There  were  several  periods  during  which  local  disturbances  generated  a  great  many  detec¬ 
tions  at  the  ARCESS  array.  These  are  very  near-field  events  that  occur  in  particular  when 
the  temperature  is  below  freezing  and  suddenly  drops  more  than  10  degrees  Celsius. 

Table  3.4.2  shows  daily  and  hourly  distribution  of  detections  for  ARCESS. 

Events  automatically  located  by  ARCESS 

During  days  091 . 1991,  through  273, 1991,  3999  local  and  regional  events  were  located  by 
ARCESS,  ba.sed  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
21 .9  events  per  processed  day  (183  days  processed).  46  %  of  these  events  are  within  300 
km,  and  86  %  of  these  events  are  within  1000  km. 

FINESA  detections 

The  number  of  detections  (phases)  reported  during  day  091  1991,  through  day  273  1991, 
was  34,925,  giving  an  average  of  275  detections  per  processed  day  (127  days  processed). 

Table  3.4.3  shows  daily  and  hourly  distribution  of  detections  for  FINESA. 
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Events  automatically  located  by  FINESA 

During  days  091  1990,  through  274  1991, 2095  local  and  regional  events  were  located  by 
FINESA,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  16.5  events  per  processed  day  (127  days  processed).  64  %  of  these  events  are  within 
300  km,  and  86  %  of  these  events  are  within  1000  km. 


GERESS  detections 

The  number  of  detections  (phases)  reported  during  day  091  1991,  through  day  274  1991, 
was  34,007,  giving  an  average  of  199  detections  per  processed  day  (171  days  processed). 

Table  3.4.4  shows  daily  and  hourly  distribution  of  detections  for  GERESS. 

Events  automatically  located  by  GERESS 

During  days  091  1991,  through  273  1991,  2783  local  and  regional  events  were  located  by 
GERESS,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  16.3  events  per  processed  day  (171  days  processed).  71  %  of  these  events  are  within 
300  km.  and  89  %  of  these  events  are  within  1000  km. 


Poland  detections 

The  number  of  detections  (phases)  reported  by  the  station  KSP  during  day  176  1991, 
through  day  273  1991,  was  38,143,  giving  an  average  of  397  detections  per  processed  day 
(%  days  processed). 

The  data  transmitted  from  the  station  SFP  were  in  general  very  unreliable  throughout  the 
reporting  period,  with  lots  of  gaps,  wrong  gain  and  wrong  time.  It  was  not  possible  to  pro¬ 
cess  the  data  received  from  the  SFP  station  during  this  period. 


J.  Fyen 
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3.5  IMS  operation 

The  Intelligent  Monitoring  System  (IMS)  was  installed  at  NORSAR  in  December  1989 
and  was  operated  at  NORSAR  from  1  January  1990  for  automatic  processing  of  data  from 
ARCESS  and  NORESS.  A  new  version  of  IMS  that  accepts  data  from  an  arbitrary  number 
of  arrays  and  single  3-component  stations  has  now  been  developed  by  SAIC,  and  has  been 
operated  at  the  Center  for  Seismic  Studies  (CSS)  in  Rosslyn,  VA,  since  March  1991.  This 
new  version  will  be  installed  at  NORSAR  in  October  1991 .  In  anticipation  of  the  new  ver¬ 
sion  of  IMS,  the  operation  at  NORSAR  of  the  first  IMS  version  was  discontinued  on  30 
June  1991.  During  July  1991,  NORSAR  personnel  received  a  training  course  at  the  CSS 
for  operation  of  IMS  Version  2  and  are  well  prepared  for  taking  over  the  responsibility  for 
the  daily  operation  of  IMS  in  October  this  year. 

The  statistics  in  the  following  reflect  IMS  operation  at  NORSAR  during  April  to  June 
1991.  In  general,  our  routine  operation  of  IMS  during  this  period  progressed  well,  and  the 
system  proved  to  be  vei7  powerful  and  flexible.  The  well-developed  automatic  processing 
combined  with  very  versatile  interactive  tools  has  kept  the  analyst  workload  at  a  low  level, 
and  in  fact  only  one  analyst  was  required  to  handle  the  regular  processing  at  any  lime. 

Table  3.5. 1  presents  an  overview  of  IMS  event  processor  downtimes.  Table  3.5.2  gives  a 
summary  of  phase  detections  and  proce.ssed  regional  events  by  IMS.  From  top  to  bottom, 
the  table  gives  the  total  number  of  detections  by  the  IMS,  the  number  of  detections  that  are 
associated  with  events  declared,  the  number  of  detections  that  are  not  associated  with  any 
events,  the  number  of  regional  events  declared  by  the  IMS,  the  number  of  such  events 
rejected  by  the  analyst,  the  number  of  events  accepted  by  the  analyst,  the  number  of  events 
accepted  by  the  analyst  without  any  changes,  and  finally  the  number  of  events  accepted 
after  some  sort  of  modifications  by  the  analyst.  This  last  category  is  divided  into  two 
classes:  Events  where  pha.ses  (not  delected  by  the  IMS)  have  been  added  by  the  analyst 
and  events  for  which  the  phase  assignments  by  the  IMS  have  been  changed  or  one  or  more 
phase  detections  have  been  removed. 

From  initial  review  of  the  IMS  bulletin,  it  is  clear  that  the  final  output  is  of  very  high  qual¬ 
ity  from  a  .seismological  point  of  view.  The  IMS  was  u.sed  extensively  for  the  reporting  of 
regional  events  during  the  GSETT-2  test,  conducted  by  the  Geneva  Group  of  Scientific 
Experts  during  22  April  -  2  June  1991. 


B.Kr.  Hokland 
U.  Baadshaiig 
S.  Mykkeitvcit 
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1991  Apr  14  19;3();00.000-  1991  Apr  14  24;0():(X).(X)0 
1991  Apr  17  15:30:00.000-  1991  Apr  18  02:(X):(K).000 
1991  Apr  18  17:00:00.000-  1991  Apr  19  07:00:00.000 
1991  Apr  19  19:00:00.000-  1991  Apr  20  05:00:00.000 
1991  Apr  20  18:00:00.000-  1991  Apr  22  24:00:00.000 
1991  Apr  26  17:00:00.000-  1991  Apr  27  05:30:00.000 


Table  3.5.1.  IMS  event  prex’essor  (pipeline)  downtimes. 


Apr  91 

May  91 

Jun  91 

Total 

Phase  detections 

16714 

17063 

16565 

50342 

-Associated  phases 

2849 

2656 

2392 

7897 

-Unassociated  phases 

13865 

14407 

14173 

42445 

Events  automatically  declared  by  IMS 

1080 

1183 

973 

3236 

-Rejected  events 

197 

286 

176 

659 

-Accepted  events 

883 

897 

797 

2577 

Unchanged  events 

296 

227 

229 

752 

Mtxiified  events 

587 

670 

568 

1825 

Phases  added 

67 

49 

69 

185 

Phases  changed  or  removed 

520 

621 

499 

1640 

Table  3.5.2.  IMS  phase  detections  and  event  summary. 
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4  Improvements  and  Modifications 

4.1  NORSAR 

NORSAR  data  acquisition 

No  modification  has  been  made  to  the  NORSAR  data  acquisition  system.  There  is  an 
increasing  problem  with  loss  of  data. 

NORSAR  detection  processing 

The  NORSAR  detection  processor  has  been  running  satisfactorily  on  the  IBM  during  this 
reporting  period. 

Detection  statistics  are  given  in  section  2. 


NORSAR  event  processing 

There  are  no  changes  in  the  routine  processing  of  NORSAR  events,  using  the  IBM  sys¬ 
tem. 


NORSAR  refurbishment 

The  process  of  evaluating  and  testing  technical  options  for  refurbishment  of  NORSAR  has 
continued.  A  contract  starting  1  Oct  1991  that  will  provide  funding  for  such  a  refurbish¬ 
ment  has  now  been  signed.  Technical  progress  will  be  reported  in  forthcoming  Semian¬ 
nual  Technical  Summaries. 

4.2  Regional  arrays  and  3-coniponent  stations 

Detection  processing 

The  routine  detection  processing  of  the  arrays  is  running  satisfactorily  on  each  of  the 
array’s  SUN-3/280  acquisition  systems.  The  same  program  is  used  for  NORSAR,  NOR- 
ESS,  ARCESS,  FINESA,  GERESS,  KSP,  SFP  but  with  different  recipes.  The  beam  table 
for  NORESS  and  ARCESS  is  found  in  NORSAR  Sci.  Rep.  No.  1-89/90.  The  beam  table 
for  FINESA  and  GERESS  is  found  in  NORSAR  Sci.  Rep.  No.  1-90/91. 

Detection  statistics  are  given  in  section  3. 
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Event  processing.  Phase  estimation 

This  process  will  do  F-K  and  polarization  analysis  for  each  detection  to  identify  phase 
velocity,  azimuth  and  type  of  phase,  and  the  results  are  put  into  the  ORACLE  detection 
data  base  for  use  by  IMS.  Data  for  the  three-component  Polish  stations  have  been  made 
available  for  IMS. 

Plot  and  epicenter  determination 

A  description  of  single  array  event  processing  is  found  in  NORSAR  Sci.  Rep.  No.  2-88/ 
89,  and  NORSAR  Sci.  Rep.  2-89/90. 


J.  Fyen 
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5  Maintenance  Activities 

5.1  Activities  in  the  field  and  at  the  Maintenance  Center 

This  section  summarizes  the  activities  in  the  field,  at  the  Maintenance  Center  (NMC) 
Hamar,  and  NDPC  activities  related  to  monitoring  and  control  of  NORSAR,  NORESS  and 
ARCESS.  Activities  at  other  field  sites  are  also  listed.. 

Activities  involve  preventive  and  corrective  maintenance,  modification  of  equipment,  for 
example,  in  connection  with  interfacing  to  the  new  NORSAT  B  satellite  station  (RW. 
Larsen,  ARCESS,  August  1991). 

NORSAR 

This  array  was  visited  in  May,  June,  July,  Augu.st  and  September  1991.  As  in  previous 
periods,  activities  related  to  this  array  have  been  diverse,  involving:  preventive  and  cor¬ 
rective  maintenance  such  as  line  and  communications  tests,  modem  repair,  SLEM  test, 
adjustment  of  mass  position  (MR)  and  free  period  (FP)  LP  instruments,  adjustment  of 
channel  gain  and  DC  offset  for  all  SP/LP  channels,  replacement  of  LP  seismometer 
Remote  Centering  Devices  (RCDs),  cable  splicing,  RA-5  amplifier  adjustments  and 
replacements,  replacement  of  main  fuses  after  lightning,  replacement  of  protection  cards, 
replacement  of  discrete  output  card. 

NORESS 

The  NORESS  array  was  visited  in  June  and  August  1991  in  connection  with  Hub  unit 
replacement  after  serious  damage  by  lightning,  and  repair  of  Hub  1/F  card  which  did  not 
send  synch,  pulses  to  remote  sites. 

ARCESS 

This  array  was  visited  in  April,  May.  June,  July,  August  and  September  -  in  connection 
with  activities  such  as:  UPS  failure  after  power  outage,  replacements  of  fiber  optical 
receiver/transmitter  at  sites  A2  and  A3,  modifications  of  all  seismometer  vaults,  cleaning 
of  fiber  optical  connectors,  modification  of  "white  box"  to  achieve  correct  I/F  to  the  new 
NORSAT  B  satellite  station,  adjustment  of  fiber  optical  link  to  A2  and  D6,  removal  of  the 
old  NORSAT  BO  telemetry  rack  from  the  Hub. 

FINESA 

FINESA  was  visited  in  April,  July  and  September  in  connection  with  field  maintenance, 
replacement  of  RA-5  amplifier,  and  replacement  of  cable  to  remote  sites  after  lightning. 

Poland,  Ksiaz/Stary  Foiwark 

Problems  which  had  been  present  ever  since  installation  of  the  station  were  partially 
solved  during  a  visit  in  April,  and  the  communication  link  tested. 
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Subarray/ 

Area 

Task 

Date 

Poland/Ksiaz/ 
Stary  Folwark 

Representatives  from  Teledyne  Geotech,  NTA 
satellite  department  and  NORSAR  solved  problems 
which  had  been  present  ever  since  the  installation 
of  the  stations.  During  communications  tests,  both 
NTA/Lillestram  and  NORSAR  were  engaged. 

13-20  Apr 

ARCESS 

Local  representative  restored  the  operation  after 

UPS  failure 

17  Apr 

FINESA 

NMC  (0.  Loken)  carried  out  field  maintenance 

15,16,17  Apr 

NDPC 

Daily  check  of  all  aiTays,  i.e.,  NORSAR,  NORESS, 
ARCESS,  FINESA.  GERESS  and  Poland  (KSP.SFP). 
NORSAR  SP/LP  instruments  calibrated  every  week. 
Free  Period  (FP)  and  Mass  Position  (MP)  have  been 
measured  and  adjusted  when  outside  tolerances. 

April 

NORSAR 

04C 

Communications  test.  Replaced  AHS-card  in  CTV 
modem 

23  May 

ARCESS 

Replaced  fiber  optical  receiver/transmitter  at  sites 

A2  and  A3. 

14,15,16 

May 

NDPC 

Daily  check  of  all  arrays,  i.e.,  NORSAR,  NORESS, 
ARCESS,  FINESA,  GERESS  and  Poland  (KSP.SFP). 
NORSAR  SP/LP  instruments  calibrated  every  week. 
Free  Period  (FP)  and  Mass  Position  (MP)  have  been 
measured  and  adjusted  when  outside  tolerances. 

May 

Table  5.1.  Activities  in  the  field  and  the  NORSAR  Maintenance  Center,  including  NDPC 
activities  related  to  NORSAR.  NORESS.  ARCESS,  FINESA,  GERESS  and  Poland  (KSP, 
SFP)  arrays.  1  April  -  30  Sepember  1991. 
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Subarray/ 

Area 

Task 

Date 

NORSAR 

02B 

Line  and  SLEM  tests  carried  out 

27  June 

ARCESS 

All  seismometer  vau'  :  modified 

10-21  June 

NORESS 

Hub  unit  replaced  after  serious  damage  by  lightning 

18-25  June 

Apatity, 

Installed  a  Nanometrics  data  acquisition  system 

10-15  June 

USSR 

unit  and  a  3-component  S 13  station 

NDPC 

Daily  check  of  all  arrays,  i.e.,  NORSAR,  NORESS, 
ARCESS,  FINESA,  GERESS  and  Poland  (KSP.SFP). 
NORSAR  SP/LP  instruments  calibrated  every  week. 
Free  Period  (FP)  and  Mass  Position  (MP)  have  been 
measured  and  adjusted  when  outside  tolerances. 

June 

NORSAR 

OlA 

Main  fuses  replaced  after  lightning.  Comm,  line 
quality  tested.  Adjusted  MP/FP  all  LP  instruments. 
Adjusted  channel  gain  and  DC  offset  all  SP/LP 
channels 

16  July 

GIB 

Replaced  MP  and  RP  Remote  Centering  Devices 
(RCDs)  NS  LP  seismometer.  Adjusted  all  MP/FP 

RCDs.  Adjusted  channel  gain  all  SP/LP  channels 

2  July 

02  B 

Replaced  MP/FP  RCDs  EW  seismometer.  Replaced 
discrete  output  card  (XI). 

19  July 

Adju.sted  MP/FP  RCDs  all  LP  instruments 

15  July 

02C 

Adjusted  MP/FP  RCDs  all  LP  instruments.  Adjusted 
channel  gain  and  DC  offset  all  SP/LP  channels 

18  July 

03C 

Line  check  B/C-loop 

1  July 

Table  5.1  (cont) 
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Subarray/ 

Area 

Task 

Date 

04C 

Adjusted  MP/FP  RCDs  all  LP  instruments.  Adjusted 
channel  gain  and  DC  offset  all  SP/LP  channels 

1  July 

06C 

Replaced  protection  card  SP03.  Adjusted  MP/FP 
all  LP  seismometers.  Adjusted  channel  gain  and 

DC  offset  all  SP/LP  channels 

17  July 

ARCESS 

Replaced  fiber  optical  transmitters  D4 

Channel  connectors  fiber  optical  cables  A3 

25-26  July 

FINESA 

Replaced  all  RA-5  amplifiers 

9,10,11  July 

NDPC 

Daily  check  of  all  arrays  has  been  carried  out  weeks 

27, 28  and  29.  i.e.,  NORSAR,  NORESS,  ARCESS, 
FINESA.  GERESS  and  Poland  (KSP.  SFP).  The 
actual  weeks  NORSAR  LP  instruments  were  cali¬ 
brated,  Free  Period  (FP)  and  Mass  Position  (MP) 
were  measured  and  adjusted  (those  outside  speci¬ 
fications). 

July 

NORSAR 

OIB 

Adjusted  RA-5  amplifier  SP02 

19  Aug 

02C 

Comm,  check  in  cooperation  with  NTA/Lillestrom 

19  Aug 

02B 

Remote  Centering  Devices  (RCDs)  MP/FP  all 

LP  seismometers  replaced 

21-22  Aug 

06C 

Adjusted  RA-5  amplifier  SP02.  RA-5  amplifier  SP05 
replaced. 

23  Aug 

NORESS 

Hub  I/F  card  repaired  as  it  did  not  send  synch  pulses 
to  the  remote  sites. 

2  Aug 

Lightning  caused  array  outage  and  deterioration  of 
equipment  throughout  the  entire  array,  such  as;  Hub, 
remote  sites,  TOD  (clock),  terminals  (2),  converter 
and  modems. 

16  Aug 

Table  5^  (cont) 
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Subarray/ 

Task 

Date 

Area 

Array  resumed  operation,  but  sites  A1 ,  C2,  D2  and 

D8  were  still  inoperative. 

27  Aug 

All  sites  performed  satisfactorily,  but  A1  was  still 
inoperative  due  to  lack  of  spare  parts  for  the  CPU. 

30  Aug 

ARCESS 

“White  box"  modified  during  the  stay.  Afterwards 
succes.sfully  interfaced  to  the  new  NORSAT  B 
satellite  station. 

13-16  Aug 

Svalbard 

Site  survey  and  measurements  carried  out  in  order 
to  find  a  suitable  location  for  the  planned  miniarray 
in  the  Longyearbyen  area. 

19-21  Aug 

NDPC 

Daily  check  of  all  arrays  has  been  carried  out,  i.e., 
NORSAR,  NORESS,  ARCESS,  FINESA,  GERESS, 
and  Poland  (KSP,  SFP).  Every  week  the  NORSAR 

SP/LP  instruments  were  calibrated.  Free  Period  (FP) 
Mass  Position  (MP)  were  measured  and  adjusted 
(those  outside  specifications). 

August 

NORSAR 

OIB 

Replaced  RA-5  amplifier  SP02  vault.  Cable  spliced 

SP04.  Adjusted  FP/MP  all  LP  instruments 

23,25  Sep 

ARCESS 

Adjusted  fiber  optical  links  for  sites  A2  andD6. 

Removed  all  NORSAT  BO  telemetry  racks  from  Hub 

14-16  Sep 

FINESA 

Replaced  cables  to  remote  sites  after  lightning,  for  sites 
A2,  B 1 ,  B2,  B3,  B6,  Cl ,  C4,  C5  and  C6. 

10-13  Sep 

NDPC 

Daily  check  of  all  airays  has  been  ciuried  out,  i.e., 
NORSAR.  NORESS,  ARCESS,  FINESA,  GERESS, 
and  Poland  (KSP,  SFP).  Every  week  the  NORSAR 
SP/LP  instruments  were  calibrated.  Free  Period  (FP) 
Mass  Position  (MP)  were  measured  and  adjusted 
(those  outside  specifications. 

September 

Table  5.1  (cont.) 
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5.2  Array  status 

As  of  30  September  1991  the  following  NORSAR  channels  deviated  from  tolerances: 

01 A  01  8  Hz  filter 

02  8  Hz  filter 

04  30  dB  attenuation 


O.A.  Hansen 
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6  Documentation  Developed 

Kvaema,  T.  and  DJ.  Doombos,  Scattering  of  regional  Pn  by  Moho  Topography,  Geophys. 
Res.  Lett.,  Vol.  18. 1273-1276, 1991. 

Kvaema,  T.  aind  F.  Ringdal,  Integrated  array  and  3-component  processing  using  a  seismic 
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7  Summary  of  Technical  Reports  /  Papers  Published 

7.1  Statistical  optimization  of  seismic  holography  algorithms  for 
array  data  processing 

In  this  paper  we  describe  the  application  of  statistical  parameter  estimation  theory  to  the 
problem  of  locating  weak  seismic  radiation  sources  in  the  lithosphere.  The  radiation  may 
have  two  origins.  The  first  is  scattering  generated  by  an  earthquake  or  other  type  of  source. 
These  scattered  seismic  waves  are  to  be  extracted  from  the  background  primary  waves. 
The  other  origin  is  a  weak  seismic  emission  in  the  medium.  The  problem  of  mapping  of  a 
weak  seismic  radiation  source’s  spatial  distribution  on  the  basis  of  seismic  array  data  pro¬ 
cessing  is  called  seismic  holography. 

Theory 

The  seismic  field  §  (i,  p)  measured  on  the  surface  may  be  given  as  the  sum  of  the  “signal” 
field,  produced  by  the  source,  located  in  the  medium  at  point  r,  and  the  residual  field 
n('.p)  : 


K^P) 


+  oo 


J  H5(t-T)C(T,^p)</T+T|(r,  p) 


(1) 


where  p  is  the  point  on  the  surface,  ju  (/)  is  a  scalar  waveform  generated  by  the  source  at 
point  r,  G  (t,  r,  p)  is  the  Green  function  of  the  medium,  and  p  is  a  scaling  factor.  Accord¬ 
ing  to  eq.  ( 1 ),  the  multidimensional  time  series  x,  =  p, ) . | ( i/f^, p„) )  (x,  is  a 

column  vector)  recorded  by  one-component  array  sensors  has  the  following  structure 


T  =0 

where  J,  =  h,Cr)  =  G(T,r,p,),/6  i,m  is  a  column  vector  of  the  medium’s 

impulse  response  for  the  seismic  wave  propagating  from  the  source  at  point  ^  to  the  sen¬ 
sors  at  points  p,;  ii,  is  the  “noise”  time  series,  generated  by  other  sources  located  away 
from  point  r;  is  the  sampling  rate  and  *  denotes  convolution.  In  the  case  of  scattering,  it, 
contams  strong  components  generated  by  the  primary  wave,  and  may  therefore  influence 
the  estimation  of  scattered  wave  power. 

In  most  cases  there  is  no  a  priori  information  about  the  source  waveform  and  the  “noise” 
field  features.  So  from  a  statistical  point  of  view  it  is  reasonable  to  assume  the  scalar  pro¬ 
cess  r,  and  the  vector  proce.ss  to  be  realizations  of  Gaussian  stationary  time  series  with 
zero  mean  and  power  spectral  densities  0  (X.)  and  F  (X.) ,  respectively.  Under  the  additional 
J  assumption  that  and  r|,  are  statistically  independent,  the  matrix  spectral  density  of  the 

time  series  x,  becomes; 
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F^{X,B)  =e(^{X)h,a)h,  a)+F^{X),  Xe  [0,271] 

0  =  )i^ 

2n  (3) 

0 


where 

ac> 

hr(X)  =  y  h,Cr)e'"' 

I  =  0 


is  the  vector  respemse  of  the  propagation  paths  from  point  r  to  the  amiy  sensors;  G  is  the 
power  of  the  seismic  source  at  point  r  and  /  =  2n//f^  is  a  normalized  freciuency. 

For  simple  reference  models,  /(,  (>.)  may  be  computed  by  solving  the  “direct"  seismic 
problem.  With  reference  model  of  the  medium,  we  mean  the  model  which  determines  the 
main  features  of  the  seismic  wave  propagation  without  accounting  for  the  inhomogenei¬ 
ties  which  are  to  be  detected  by  the  analysis  of  the  scattered  waves.  This  means  that  in 
etjuation  (2).  we  ignore  the  effect  of  secondary  .scattering  of  the  .seismic  waves  caused  by 
the  inhomogeneities. 

Ei|uations  (2)  and  (3)  allow  consideration  of  the  problem  of  locating  seismic  scatterers  as 
the  statistical  problem  of  source  power  0  estimation,  which  is  successively  solved  for 
each  pennt  of  the  "scanned"  medium  area  (Troitsky  e/  al,  1980).  Under  the  assumption 
that  the  power  spectra  and  F  (> )  are  known,  statistically  optimal  algorithms  may  be 
proposed  for  this  estimation  problem.  But  in  practice  these  assumptions  are  never  fulfilled 
and  these  algorithms  must  be  modified  to  make  them  efficient  also  when  we  do  not  have 
full  information  about  statistical  features  of  the  observations. 

Statistically  optimal  algorithms  for  estimation  of  the  seismic  source  power  0  may  be 
developed  using  the  maximum  likelihood  method  (in  its  asymptotical  modification), 
llsing  the  Bartlett  formula  (Bartlett,  1951 )  for  inversior  of  the  matrix  in  equation  (3),  we 
obtain: 


F 


F:'/;/7  f 


-t 


0(J) 


1  +h 


F"'/) 


(4) 


For  the  maximum  likelihood  estimate  (\  of  the  power  of  a  seismic  source,  located  at  point 
f,  Kushnir  ( 1989)  shows  that  the  following  eijuation  holds: 


.V 


I 

I  -  1 


0  V'‘ 

V  ./ 


+  0vO  V  I 

'  I  /  >-- 
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(5) 
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where  =  h,' 0.j)Xj,  =  Vik^),  V(k)  =  hr^('k)F^^  ('k)hr(k) , 

s 

Xj  =  ( 1  /  (  Jn)  )  ^  .  }  =  2nj/N,  N  is  the  number  of  observations,  and  t  denotes  the 

j=  I 

Hermitian  conjugation. 

Equation  (5)  must  in  general  be  solved  by  numerical  methods,  but  in  some  practically 
imporutnt  cases  the  solution  may  be  written  in  “closed”  form.  Namely,  if  the  source 
power  is  much  less  than  the  "noise”  power,  then  from  (5)  the  following  approximate  solu¬ 
tion  may  be  obtained  (Kushnir,  1989): 


j=  1 


This  is  a  consistent  and  asymptotically  normal  estimate,  but  it  has  larger  variance  than  the 
exact  solution  of  et|uation  (5). 

listimates  based  on  solving  eci'aation  (.‘i)  do  not  provide  answers  to  the  problem  of  map¬ 
ping  of  scattered  seismic  radiation  sources  formulated  in  the  introduction,  because  such 
C'-iimates  assume  that  the  povser  spectnim  F„t/ )  of  the  “noise"  and  source  waveform 
povu'r  spectrum  are  known.  It  can  be  shown  that  the  information  about  is 

\er\  import.int  because  F,,  ( /  '  determines  the  statistical  propenies  of  all  power  H  estima¬ 
tors.  W  hen  K„  /  /  i  is  unknown,  it  is  reasonable  t()  use  the  adaptive  approach.  In  that  case. 
F.,!/  I  in  equations  (5)  .md  (6)  is  substituted  by  F^O  )  derived  from  the  obseI^'ations 
themselves.  Hut  in  the  lirst  of  the  problems  outlined  in  the  introduction,  it  is  impossible  to 
consistently  estimate  F,,  ( /  i  from  obseivaiions.  Often  the  same  is  true  for  the  second 
problem,  .Substituting  for  F,,  (/  )  in  the  expressions  (5)  and  (6).  any  consistent  estimate 
F  1/  )  of  the  m.itrix  powei  sjiectrum  of  the  observations  1,  yields,  values  and 
.tpproachmg  /ero  m  prob.ihilitv.  In  other  words,  adaptation  of  the  estimates  in  expressions 
<  5 )  .md  (6)  IS  impossible. 

Nevcitheless.  the  solution  of  the  problem  under  the  assumption  of  full  a  priori  information 
about  the  obserx  ations'  statistic.d  features  is  useful.  It  allows  us  to  find  the  important  sta¬ 
tistic.  which  is  contained  in  all  statistically  optimal  algorithms  for  anay  data  processing. 
This  statistic  is  the  output  of  the  optimal  group  filter  (OGF): 

y,  =  'r*^,  Yj  =  R{Xj)Xj 

where 
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Raj) 


A^a.)F-'(x.) 

hUxj)F-Uij)haj) 


2n 

r,=  (1/(2ji)) 

0 


(7) 


If  the  input  time  series  x,  fits  the  observational  model  (2),  (3),  then  it  is  easy  to  show  that 
the  filter  R  (X.)  minimizes  the  mean  square  value  of  the  output  noise  component  provided 
the  signal  component  is  undistorted.  The  remarkable  feature  of  the  filter  R  (K)  is  that 
replacement  in  the  expression  (7)  of  the  unknown  spectral  matrix  (X)  of  the  noise  by 
the  spectral  matrix  F,(X)  of  the  observations,  does  not  change  its  output:  it  keeps  the  sig¬ 
nal  undistorted  and  minimizes  the  noise  component  power.  The  new  filter 


KiX)  -  - - 

h  '(X)f;’  a)  ha) 


(8) 


is  equivalent  to  the  filter  given  by  equation  (7).  This  can  be  shown  using  Bartlett's  formula 
(4).  Due  to  this  property,  each  statistic,  depending  only  on  the  output  of  OGF,  has  the 
adaptive  feature.  This  means  that  in  this  stati.stic  we  may  substitute  instead  of  the  unknown 
noise  matrix  power  spectrum  F^(X,)  any  consistent  estimate  F;if(X)  of  the  matrix  power 
spectrum  of  the  observations  .1,,  thus  providing  the  adaptive  (AOGF)  filter  k  (X) . 

Let  us  now  consider  some  practically  interesting  modifications  of  OGF  K  (X) .  In  the  case 
when  the  noise  field  mainly  consists  of  coherent  primary  waves,  generated  by  a  seismic 
source,  the  matrix  spectium  F^  (X)  may  be  written  in  the  form: 


F^(X)  =  eF^ (X) +\j;(X)4 (^)  4  (9) 

where  F^  (X)  is  the  matrix  power  spectrum  of  the  “diffuse  noise"  component,  e  is  a  scal¬ 
ing  factor,  v(X)  is  the  scalar  power  spectrum  of  the  primary  waveform,  and 
^  (X)  =  (<?,(X),l  e  1.  m)  is  the  frequency  response  of  the  medium  for  the  primary  wave 
paths  from  the  source  to  the  an  ay  sensors.  Using  the  Bartlett  formula  (4)  it  is  easy  to  show 
that  for  the  matrix  spectrtim  (9),  OGF  K(X)  has  the  form: 


R{k)  =  R„a)  +Oy(Z), 


Roa) 


P  (X)B(X) 
(X)B(X)MX) 


(10) 


where 


B(X) 


V  1^(X)2|  / 


/  is  the  identity  matrix. 


(7x.(e)||  ->0  if  e-xO 
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!!)"(>.)  II  is  the  norm  of  the  vector  function 7.  So  if  the  “diffuse”  components  of  the  noise 
and  signal  fields  are  weak  and  hence  the  matrices  F“‘  (k)  and  F”‘  (k)  are  close  to  singular, 
it  is  reasonable  to  use  the  stable  form  of  OGF  with  frequency  response  Ro(k) . 

In  the  opposite  particular  case,  when  the  noise  field  is  purely  “diffuse”,  i.e.,  it  may  be 
modeled  as  white  noise,  the  matrix  power  spectrum  F^(^)  =  o^(X.)I  and 

RM  =  =  f  111) 

|/>(X)| 

The  undistorted  group  filter  r(>.)  is  the  mathematical  expression  for  the  conventional 
algorithm  of  seismic  holography  (Troitsky,  et  al,  1980),  comprising  the  “focusing  of  rays” 
radiated  from  a  given  point  r  of  the  scanned  medium  area.  It  compensates  for  the  phase 
and  amplitude  distortions  of  the  airay  records  caused  by  the  wave  propagation  from  the 
source  to  the  array  sensors. 

Let  us  consider  the  dispersion  ^ ^  V  v?  of  the  group  filter  K,  (/.)  output  y,.  i  e  T7^.  In 

^1=1 

case  the  input  v,  has  a  matrix  power  spectrum  described  by  the  formulas  (3)  and  (9).  the 
property  of  the  matrix  B(  / ):  B(>  )</(>,)  £3 0  yields: 

P-  lim  di  =  0  +  P.  (12) 

u  here  -*  0  if  s  ^  0.  So  if  the  "diffuse"  noise  (caused,  particularly,  by  scatterers  located 
away  from  point  )■)  is  small  enough,  the  statistic  S;  provides  a  consistent  estimate  of  the 
scattered  waive  source  power,  thus  suppressing  the  array  noise  generated  by  the  primaiy 
wave.  The  OGF  Ki/ )  (8)  has  the  same  important  feature.  And,  finally,  it  is  clear  that  the 
relation  (  12)  is  valid  for  the  adaptive  filter  A'O  ). 

Thus,  m  the  case  of  coherent  noise,  the  problem  of  consistent  estimation  of  the  seismic 
wave  power  radiated  from  point  r  can  be  solved  rather  effectively  by  using  OGF  or 

AOfJFA'G), 

•At  the  same  time,  the  conventional  holographic  group  filter  f(>.)  cannot  suppress  the 
coherent  noise  and  does  not  provide  a  consistent  estimator  of  the  seismic  wave  power 
radiated  from  point  r.This  is  so  because  the  dispersion  of  the  holographic  filter  f  (A)  out¬ 
put  (when  the  input  v,  matrix  power  spectrum  is  defined  by  (3)  and  (9))  does  not  tend  to 
the  seismic  wave  power  radiated  from  the  point  r  when  the  “diffu.se”  noise  component 
decreases. 

Results  of  simulation 

The  benefits  of  using  the  optimal  statistical  algorithms  for  mapping  of  scattered  wave 
sources  have  been  analyzed  by  computer  modeling.  The  records  of  a  Hindu  Kush  earth¬ 
quake  on  22  February  1972  with  magnitude  5.6  recorded  by  NORSAR  have  been  used  in 
this  modeling.  To  check  the  algorithm's  capability  to  process  data  recorded  at  arrays  with 
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a  small  number  of  sensors  and  a  small  aperture,  the  six  sensors  of  the  central  NORSAR 
subarray  01 A  have  been  chosen.  The  location  of  these  sensors  are  shown  in  Fig.  7.1.1,  and 
a  record  of  the  Hindu  Kush  earthquake  is  shown  in  Fig.  7.1.2. 

The  numerical  experiments  consisted  of  processing  data  from  the  six  sensors  of  the  central 
NORSAR  subarray  01  A.  The  data  were  composed  as  a  sum  of  primary  and  scattered 
waves.  As  primary  waves  we  have  tested  two  cases:  1 )  observed  array  records  of  the 
Hindu  Kush  earthquake;  2)  A  plane  P-wave  with  an  apparent  velocity  and  azimuth  corre¬ 
sponding  to  that  of  the  Hindu  Kush  earthquake.  The  modelling  of  this  plane  wave  was 
done  to  check  how  the  deviation  from  a  plane  wavefront  affects  the  quality  of  the  seismic 
holography  algorithms.  The  scattered  wave  was  simulated  as  a  P-wave  radiated  from  a 
point  source  located  in  a  homogeneous  medium  4  km  beneath  the  NORSAR  central  subar¬ 
ray.  The  source  waveform  i,  of  the  simulated  scattered  wave  was  synthesized  as  a  sample 
function  of  the  stationary  Gaussian  random  process  with  power  spectrum  estimated  using 
the  Hindu  Kush  earthquake  P-wave  record. 

In  Figs.  7. 1.3-7. 1.5  the  results  from  processing  four  different  data  sets  are  depicted.  The 
array  data  .sets  were  the  following: 

a.  The  simulated  plane  wave  from  the  Hindu  Kush  direction 

b.  The  sum  of  a)  and  the  simulated  .scattered  wave  from  the  point  source  beneath  NOR¬ 
SAR 

c.  The  original  Hindu  Kush  earthquake  records 

d.  The  sum  of  c)  and  the  simulated  scattered  wave  fro;n  the  point  source  beneath  NOR¬ 
SAR. 

The  power  of  the  scattered  wave  was  equal  to  that  of  the  primary  wave. 

In  each  of  the  four  experiments  the  array  signals  were  processed  by  the  following  algo¬ 
rithms:  a)  the  conventional  algorithm  of  seismic  holography  based  on  the  filter  f(k)  (11); 
b)  the  statistical  algorithm  based  on  the  adaptive  group  optimal  filter  k (A.)  (8)  with  the 
power  spectral  density  inverse  matrix  Fj‘  (A.)  e.stimated  by  the  multidimensional  autore¬ 
gressive  modeling  of  the  observations  ,v, ;  c)  the  statistical  algorithm  based  on  the  undis¬ 
torted  group  filter  RgO.)  (10),  in  which  the  vector  q(k)  was  chosen  in  correspondence 
with  the  plane  wave  with  the  azimuth  and  apparent  velocity  of  the  Hindu  Kush  earthquake 
P-wave. 

The  data  processing  procedure  consisted  of  scanning,  with  steps  of  0.8  km  an  8  x  8  km 
plane  parallel  to  the  surface  containing  the  scattered  wave  source  in  the  central  point.  For 
each  point  scanned,  the  medium  frequency  response  vector  hr  (A.)  was  computed  under  the 
assumption  of  medium  homogeneity.  By  this,  the  adjustment  of  the  group  filters  for  the 
extraction  of  the  scattered  wave  radiated  from  the  given  point  was  provided.  Then  the  fil¬ 
ter  output  signals  and  their  dispersion  estimates  were  computed.  As  was  pointed  out,  the 
dispersion  of  the  outputs  of  all  filters  under  investigation  may  be  used  as  the  estimates  of 
the  scattered  wave  power  radiated  from  the  given  point.  These  estimates  are  biased 
because  of  primary  wave  power  "leakage",  and  these  biases  have  different  values  for  the 
different  filters. 
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The  number  of  samples  of  the  input  vector  time  series  .1,  was  as  a  rule  equal  to  128  in 
these  experiments,  and  conesponded  to  the  Hindu  Kush  eiuthquake  P-wave  record  length. 
The  matrix  power  spectrum  F,(X)  of  the  observations  .v,  used  for  the  undistorting  group 
filter  adaptation  was  estimated  by  multivariate  AR-modeling  of  order  5. 

In  Figs.  7.1. 3-7.1 .5  the  scattered  radiation  power  maps  are  depicted.  They  are  calculated 
as  the  results  of  the  numerical  experiments  described,  with  the  conventional  holographic 
filter  (11)  (Fig.  7.1.3),  with  the  adaptive  optimal  group  filter  (7)  (Fig.  7. 1 .4),  and  with  the 
"spatial  rejecting”  filter  (10)  which  u.ses  the  a  priori  information  concerning  the  propaga¬ 
tion  direction  of  the  primary  wave  (Fig.  7.1.5). 

In  these  figures,  the  maps  with  label  (a)  are  the  results  for  the  single  simulated  primary 
wave  processing,  with  label  (b)  the  results  for  the  sum  of  the  simulated  primary  wave  and 
the  simulated  scattered  wave  processing,  with  label  (c)  the  results  for  the  Hindu  Kush 
earthquake  P-wave  processing,  and  with  label  (d),  the  results  for  the  sum  of  the  Hindu 
Kush  earthquake  P-wave  and  simulated  scattered  wave  processing. 

Comparison  of  the  maps  in  Fig.  7. 1 .3  shows  that  the  conventional  algorithm  of  .seismic 
holography  applied  to  processing  of  data  from  anays  with  a  small  aperture  and  small  num¬ 
ber  of  sensors  does  not  provide  detection  of  the  scattered  wave.  This  may  be  explained  by 
the  effect  of  the  primaiy  wave  power  “leakage",  practically  evenly  strong  for  all  points  of 
the  scanned  area.  At  the  same  time,  as  it  is  seen  from  Figs.  7.1 .4  and  7.1.5,  the  statistical 
seismic  holography  algorithms  provide  strong  ptuiks  in  the  power  maps,  the  maxima  of 
which  coincide  w  ith  the  point  of  the  scattered  wave  source  location.  The  explanation  may 
be  that  the  mentioned  algorithms  are  capable  of  suppressing  "leakage”  by  the  mutual 
compensation  of  primary  wave  components  at  the  different  an-ay  sensors.  The  adaptive 
and  "a  priori"  algorithms  demonstrate  the  same  ability  for  the  suppression  of  the  primary 
w  ave  "leakage”  and  location  of  the  scattered  wave  source,  while  in  the  ca.se  of  detection 
of  the  scattered  wave  on  the  background  of  the  Hindu  Kush  earthquake,  the  adaptive  opti¬ 
mal  group  filter  (Fig.  7.1 .4d)  has  appemed  more  effective  than  the  "a  priori”  filter  (Fig. 

7. 1 .5d).  The  advantage  of  the  optimal  group  filter  over  an  "a  priori"  filter  in  this  experi¬ 
ment  may  be  explained  by  the  deviations  of  the  Hindu  Kush  earthquake  signal  from  a 
plane  wave.  This  deviation  significantly  reduces  the  performance  of  the  filter  «o(>T  (10). 
At  the  same  time,  for  the  adaptive  filter  A'O  )  no  assumptions  are  needed  about  the  pri¬ 
mary  wave  front  shape,  which  makes  it  more  flexible  and  effective. 
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Fig.  7.1.3.  Array  data  pnKcssing  using  conventional  holographic  group  filter,  (a)  Simu¬ 
lated  plane  wave  from  earthquake;  (b)  simulated  earthquake  plane  wave  +  simulated  scat¬ 
tered  wave;  (c)  P-wave  of  the  Hindu  Kush  earthquake;  (d)  Hindu  Kush  earthquake  P-wave 
+  simulated  scattered  wave.  The  contour  levels  reflect  the  estimated  power  of  waves  arriv¬ 
ing  from  different  points  of  the  scanned  plane.  The  power  estimates  are  not  normalized, 
hence  only  relative  values  within  each  map  should  be  considered. 
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7.2  A  new  ML  scheme  for  3-component  slowness  estimation  which 
incorporates  the  crustal  transfer  function 

Background 

Efforts  to  exploit  the  wavefield  information  contained  in  the  records  of  3-component  (3C) 
seismometers,  say  in  terms  of  polarization  characteristics,  slowness  vectors  and  other 
types  of  signal  attributes,  have  until  recently  been  less  successful  in  the  high  frequency 
band  1-100  Hz.  This  problem  has  to  a  large  extent  been  overcome  by  introducing  flexible 
particle  motion  models  in  combination  with  effective  maximum  likelihood  (ML)  schemes 
for  signal  parameter  estimation  from  3C  recordings  (Christoffersson  et  al,  1988;  Roberts 
et  al,  1989;  and  Roberts  and  Christoffersson,  1990).  Their  techniques  allow  us  to  produce 
quite  precise  automatic  preliminary  bulletins  for  single  3C  (digital)  stations  (Ruud  and 
Husebye,  1991).  A  potential  flaw  in  the  above  and  similar  techniques  is  that  the  response 
of  the  crust  beneath  the  receiver  is  not  incorporated  in  the  signal  model.  However,  recently 
Kennett  (1991)  suggested  a  new  approach  for  decomposing  the  3C  seismograms  of  R  SV 
and  SH  waves  incident  onto  the  free  surface  by  assuming  that  the  slowness  vector  and  the 
layered  crustal  structure  beneath  the  station  are  known.  A  similar  method  can  in  principle 
be  used  for  solving  the  inverse  problem,  and  this  is  the  subject  dealt  with  here. 

Method 

The  problem  is  that  of  estimating  the  azimuth  p  and  horizontal  slowness  p  of  the  plane  P- 
wave  incident  onto  the  layered  structure  below  the  free  surface.  In  mathematical  form,  the 
model  of  our  3C  recordings  can  be  written  as: 

M(a))  =  f{{Oj)h(Wj,s,m)  (1) 

where  is  the  vector  of  discrete  Fourier  amplitudes  at  angular  frequency 

(IV  for  the  individual  recording  components  (NS,  EW,  Z)  calculated  for  a  time  window 
containing  the  signal;  /(ov)  is  the  product  of  the  incoming  signal  amplitude  and  the  seis¬ 
mometer  frequency  response;  *  (o>^,  i,  m)  =  t  ^  is  the  crustal  transfer  function  for  an 

incoming  P-wave  at  frequency  /(ov)  for  slowness  vector  i (p.p) ,  and  the  presumed  known 
crustal  structure  given  by  the  model  vector  m;  a  {la.)  =  n,,  n^)  ^  is  the  noise  vector  and  T 

means  transposed.  Assuming  the  noise  is  Gaussian  and  uncorrelated  at  different  frequen¬ 
cies,  then  FAn^n-'  )  =  where  E  is  mathematical  expectation,  the  asterisk  *  implies 
Hermitian  conjugation  and  6^^  is  the  Kronecker  delta  symbol.  With  the  above  assumption 
it  can  be  shown  that  the  maximum  likelihood  (ML)  functional  will  be 

L  =  1/25;  ln(detC,.)  -  1/2  J  (2) 

J  j 

where  the  summation  is  over  N  di-screte  frequencies  in  the  selected  signal  bandwidth.  L 
depends  not  only  on  the  u.seful  (informational)  unknown  parameters  p  and  p  of  the  slow¬ 
ness  vector,  but  also  on  the  unknown  regressor  parameters  /(«.) .  The  ML  estimates 
derived  by  maximization  of  the  L  functional  (2)  with  respect  to  all  unknown  parameters 
will  for  a  constant  number  of  nuisance  parameters  /(ov)  but  an  increasing  sample  size  be 
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asymptotically  optimal  estimates  of  the  useful  slowness  vector  parameters  (Kushnir  and 
Lokshtanov,  1988). 

Presuming  i(<p.p)  is  fixed,  m  is  known  a  priori  and  that  the  noise  is  orthogonal  between 
components  and  with  equal  power  at  all  components  and  frequencies,  i.e.,  Cj  =  o^i  (/  is  the 
identity  matrix),  the  maximum  of  L  with  respect  to  fj  is  achieved  at  the  points 

fj  =hj  Uj/it*  hj  (3) 

By  substituting  the  fj  here  into  (2)  we  have  that  the  ML  estimates  of  the  informational 
parameters  can  be  obtained  by  minimizing  the  following  functional 

Q(<p,p) 

j 

Now  the  functional  (4)  depends  only  on  the  observations  i.  and  the  informational  parame¬ 
ters  <p,p. 

_  -  _ *  -  ^/2 

By  using  normalized  transfer  functions  in  the  formula  <'4),  i.e.,  gj  =  hj/  {hj  hj)  we 
obtain 


©  (<P.P)  =  X  ^ 

j 

Thus  minimizing  (4)  is  equivalent  to  maximizing  the  sum  of  squai  es  of  projections  of  our 
observational  vectors  Uj  onto  the  theoretical  directions  determined  by  the  model  vectors  gj. 

The  ratio  of  the  projected  to  the  observed  energy  in  the  selected  frequency  band 

(2L  =  J  [  (g*  iij)*  {g*  Uj)  ]  /  J  {u*  u)  (6) 

j  j 

can  be  used  as  a  quality  measure  of  our  estimation.  For  example,  if  we  have  pure  signal, 
n.  =  0  for  all  j  in  eq.  ( 1 ),  then  QL  will  be  equal  to  1 .  When  a.  increases,  QL  will  decrease. 

The  above  slowness  estimation  technique  is  described  in  detail  in  a  paper  by  Lokshtanov 
et  al  ( 1991 ),  which  also  includes  its  extension  to  an  array  of  M  3C  stations  presuming  s. 
and  n,  to  be  identical. 

Results  --  velocity  analysis  ofNORESS  events 

For  the  frequency  band  of  interest  for  our  analysis  (1-15  Hz)  the  most  important  part  of  the 
crust  would  be  the  upper  5  km  as  conversion  and/or  reflections  from  deeper  layers  will  be 
separated  in  time  from  the  primary  phase.  To  our  knowledge  the  only  available  informa¬ 
tion  on  vertical  velocity  structures  in  this  part  of  the  crust  comes  from  a  study  of  Loksh¬ 
tanov  et  al  ( 1991 ).  where  upper  crust  velocity  models  were  found  by  inversion  of  short 
period  Rayleigh  wave  fundamental  mode  phase  velocity  measurements  from  NORSAR 
subarrays.  Their  results  suggest  an  upper  low  velocity  layer  about  1  km  thick  or  alterna¬ 
tively  a  two  layer  model  over  a  half  space.  Model  parameters  are  given  in  Table  7.2.1 
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together  with  those  for  uniform  velocity  halfspace  model.  For  comparison  all  three  models 
have  been  used  in  the  analysis  of  the  event  examples,  and  in  the  result  presentation  below 
they  are  referred  to  as  the  halfspace,  one-layer  and  two-layer  models. 

Analysis  procedure:  The  analysis  was  performed  both  for  the  central  single  3C  station  and 
for  the  array  of  four  3C  stations.  After  azimuth  was  determined,  the  maximum  of  the  qual¬ 
ity  parameter  QL  and  the  corresponding  apparent  velocity  was  found  by  searching  over  a 
predefined  slowness  interval.  The  results  are  presented  for  a  sliding  time  window  with 
20%  cosine  tapering  and  with  an  updating  interval  of  1/10  of  the  window  length.  The  ref¬ 
erence  time  for  the  estimated  signal  attributes  is  for  the  center  of  the  time  window. 

Event  1 "  Teleseismic  earthquake  (Hindu  Kush).  This  event,  shown  in  Fig.  7.2. 1 ,  has  a  very 
high  SNR  and  is  dominated  by  relatively  low  frequency  energy.  The  signal  has  an  emer¬ 
gent  onset  and  it  starts  at  9.2  s  on  the  given  time  scale.  The  length  of  the  sliding  time  win¬ 
dow  was  2  s  and  the  analyzed  frequencies  are  from  0.5  to  5.0  Hz.  The  QL-values  rise  as 
soon  as  the  window  moves  into  the  signal  and  stay  close  to  1.0  well  into  the  signal.  Only 
for  the  velocity  estimation  do  we  find  significant  differences  for  the  three  models.  The 
estimates  for  the  layered  models  are  quite  stable  and  close  to  the  expected  velocity  for  an 
event  at  this  distance  ( 1 3.9  km/s).  In  contrast,  the  halfspace  models  give  estimates  biased 
towards  high  velocity  -  static  corrections  can  offset  part  of  the  bias.  Regarding  azimuth 
estimation,  both  the  single  station  and  the  array  estimates  are  biased  relative  to  the 
expected  azimuth  of  96‘’  (about  10®  for  the  array  and  15®  for  single  station). 

Event  2  -  Local  explosion  (SE  Norway).  The  window  length  is  here  1  s  and  the  analyzed 
frequency  band  is  from  6  to  12  Hz.  As  seen  from  Fig.  7.2.2  the  QL-values  are  also  much 
lower  for  this  event  and  in  particular  for  the  array  where  the  signal  correlation  between 
stations  is  almost  lost  soon  after  the  onset.  For  the  three  models  a  clear  model  dependency 
is  seen  in  the  velocity  estimation  for  the  single  station.  For  the  array  it  seems  that  the  kine¬ 
matic  properties  of  the  wavefield,  i.e.,  phase  shift  betw  :en  identical  components  at  differ¬ 
ent  stations,  serves  to  stabilize  the  velocity  estimates.  This  effect  is  hardly  seen  for  the 
teleseismic  event  where  the  low  signal  frequencies  result  in  much  smaller  phase  shifts 
between  stations.  For  the  later  part  of  the  signal  (after  10.0  s).  the  velocities  are  very  high 
both  for  the  array  and  the  single  siation.  By  vespagram  analysis  of  all  the  vertical  NOR- 
ESS  sensors,  we  also  find  high  velocities  in  this  time  interval  (9-10  km/s)  although  not  as 
high  as  for  the  3C  processing. 

Discussion 

For  the  uniform  halfspace  model  the  transfer  function  is  frequency  dependent  and  our  esti¬ 
mation  method  becomes  similar  to  previous  techniques,  e.g.,  Esmersoy  et  al  (1986). 
Although  the  inclusion  of  layered  structures  clearly  improved  the  estimated  velocities,  we 
do  not  obtain  a  corresponding  improvement  of  the  quality  parameter  QL  which  measures 
how  well  the  observed  data  fit  our  signal  model.  The  reason  for  this  seems  to  be  that  the 
layered  models  affect  mainly  the  ratio  of  the  radial  to  the  vertical  amplitudes,  while  phase 
differences  between  components  remain  small  -  everywhere  within  30  from  that  of  the 
halfspace.  From  the  given  examples  and  from  general  experience  with  the  method  it  is 
seen  that  analysis  of  high  frequency  signals  is  relatively  less  successful  than  for  low  fre- 
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quencies.  An  explanation  here  might  be  that  the  models  used  are  not  adequate  at  high  fre¬ 
quencies  —  being  derived  from  Rg  waves  at  periods  from  0.7  to  1.7  s  —  and/or  that  lateral 
variations  are  highly  significant. 

In  conclusion,  we  have  here  demonstrated  the  importance  of  incorporating  the  crustal 
transfer  in  3C  analysis  methods  to  ensure  accurate  estimates  of  P-wave  slowness  -  the 
technique  will  be  extended  to  analysis  of  SV  and  SH  waves. 

E.S.  Husebye 

B.O.  Ruud,  Oslo  University 
D.E.  Lokshtanov,  Norsk  Hydro 
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Model 

P-vel. 

(kni/s) 

S-vel. 

(km/s) 

Density 

(g/cm^) 

Thickness 

(km) 

! 

!  layer 

5.10 

2.95 

2.60 

0.92 

halfspacc 

6.14 

3.55 

2.70 

oo 

2 

2-laycr 

4.93 

2.85 

2.60 

0.5 

5.40 

3.12 

2.60 

0.5 

halfspacc 

6.12 

3.54 

2.70 

CO 

3 

halfspacc 

5.9 

3.14 

2.60 

oo 

Table  7.2.1.  The  crustal  models  used  in  analysis.  Model  1  and  2  are  based  on  Rayleigh 
1  Rg)  dispersion  analysis  as  reported  by  Lokshtanov  et  al  (1991 ). 
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Fig.  7.2.1.  Analysis  results  for  Event  1  (Earthquake,  Hindu  Kush,  July  29,  1985,  origin 
time  07.54,44..5,  backazimuth  95.7,  distance  44.2).  Processing  parameters  are  given  in 
text.  From  top  and  down  the  figure  shows  3C  seismograms  (Z,E,N),  quality  parameter, 
azimuth,  apparent  velocity  for  the  array  of  four  3C  stations  and  apparent  velocity  for  the 
single  central  3C  station.  For  quality  and  velocity  the  results  for  the  three  models  are  given 
with  different  symbols:  circle  (  O  )  for  halfspace,  triangle  (A)  for  one-layer  and  square 
(□  )  for  the  two-layer  model.  The  model  independent  azimuth  estimates  are  shown  with 
plus  (+)  for  single  station  and  diamond  (  0  )  for  the  3C  array.  For  quality  the  results  for 
single  3C  and  3C  array  are  shown  in  the  same  diagram  and  are  identified  with  “3C"  and 
“A"  respectively. 
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8  10  12 


Fig.  7,2.2.  Analysis  results  for  Event  2  (Local  explosion  SE  Norway,  July  19,  1989,  origin 
time  approx.  1 1. 17.50).  Processing  parameters  are  given  in  text.  Figure  layout  and  sym- 
bols  as  for  Fig.  7.2  J . 
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Fig.  7.2.3.  Computed  radial-to-vertical  component  amplitude  ratio  as  a  function  of  fre¬ 
quency  for  the  one-layer  over  halfspace  model  in  Table  7.2.1 .  Counted  from  the  top  of  the 
six  curves  are  for  apparent  velocities  of  6, 7,  8  10, 12  and  14  km/s  respectively.  The  corre¬ 
sponding  phase  shifts  between  compionents  are  within  30  from  that  of  a  halfspace  model 
for  all  frequencies.The  strong  frequency  dependence  of  the  amplitude  ratio  illustrates  the 
need  for  correcting  velocity  estimates  for  the  effect  of  upper  crustal  layering. 
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7.3  On-line  detection  using  adaptive  statistically  optimal  algorithms 

Background 

In  a  series  of  reports  (Pisarenko  et  al,  1987;  Kushnir  et  al,  1989, 1990, 1991),  the  use  of 
adaptive  optimal  group  filtering  (AOGF)  for  the  wideband  estimation  of  seismic  signals 
has  been  investigated.  The  method  includes  autoregressive  adaptation  to  the  current  noise 
matrix  power  spectrum  and  yields  improved  noise  suppression  compared  to  conventional 
beamforming,  especially  at  low  frequencies  where  noise  power  and  coherency  is  highest. 

The  processes  have  been  implemented  in  the  NORSAR  Event  Processor  (EP)  program 
Package  (Fyen,  1989;  Kushnir  et  al,  1991),  and  comprise  adaptive  group  filtering 
(AOGF),  optimal  detection  (OD)  and  optimal  onset  time  estimation  (OE). 

The  main  advantage  of  this  signal  estimation  technique  is  that  it  retains  an  undisturbed 
signal  waveform  over  a  wide  frequency  band  (0.2-5.0  Hz),  allowing  us  to  detect  and  iden¬ 
tify  signals  at  frequencies  where  the  noise  coherency  and  power  are  the  greatest.  More¬ 
over,  when  the  final  decision  is  to  be  made  on  whether  the  estimated  signal  is  from  an 
earthquake  or  explosion  source,  it  is  important  to  have  a  wide-band  undisturbed  estimate 
of  the  signals  for  discrimination. 

The  conventional  beamforming  as  an  estimator  of  the  signal  is  together  with  bandpass  fil¬ 
tering  a  simple  and  quick  process  which  allows  us  to  detect  most  interesting  signals.  How¬ 
ever,  due  to  high  noise  coherence  for  small  arrays,  this  estimate  of  the  signal  will  be 
distorted. 

The  scope  of  this  report  is  to  analyze  the  AOGF  and  OD  for  the  purpose  of  using  the  meth¬ 
ods  in  online  processing,  i,e.,  continuous  processing  of  all  incoming  data  for  a  small  aper¬ 
ture  array  like  NORESS  and  ARCESS. 

To  use  AOGF  we  must  do  autoregressive  adaptation  to  the  noise  using  a  certain  time  win¬ 
dow  and  a  given  apparent  velocity  and  azimuth. 

We  therefore  need  to  investigate  how  often  we  need  to  adapt  to  the  noise  and  determine 
how  many  different  aiming  points  in  slowness  space  we  need  to  be  able  to  detect  the  sig¬ 
nals.  Moreover,  since  the  optimal  detector  (OD)  is  very  sensitive,  the  threshold  to  use  for 
acceptable  false  alarm  rates  must  be  carefully  considered. 

Automatic  Regional  Array  Processing 

The  current  on-line  detection  system  uses  the  STA/LTA  detector  on  conventional  bandpass 
filtered  array  beams.  A  laige  number  of  different  slownesses  are  used,  combined  with  a 
number  of  bandpass  filters.  The  airay  beams  are  formed  using  combinations  of  sub-geom¬ 
etries  of  the  array. 
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Considerable  gain  in  noise  suppression  is  obtained  by  utilizing  the  negative  correlation  of 
the  noise  wavefield  at  cenain  frequencies  and  inter-sensor  distances  (Mykkeltveit  et  al, 
1990). 

Although  the  conventional  beamforming,  filtering  and  STA/LTA  detector  scheme  work 
very  well,  requiring  small  computer  power,  there  are  a  few  problems  that  make  it  reason¬ 
able  to  evaluate  the  more  costly  statistical  optimal  algorithms  (SOA)  as  an  alternative 
online  process  tool. 

1 .  The  STA/LTA  detector  thresholds  and  bandpass  filters  are  fixed,  so  contrary  to  SOA, 
it  may  not  take  into  account  the  temporary  variations  in  noise  coherency  and  spec¬ 
trum. 

2.  As  shown  by  Fyen  (1986)  and  Mykkeltveit  et  al  (1990),  the  conventional  beamform- 
ing  provides  less  than  Js  noise  suppression  for  frequencies  below  1  Hz.  Therefore 
the  conventional  beam  detector  is  used  for  frequencies  only  above  1  Hz. 

3.  The  statistical  assumptions  for  the  STA/LTA  detector  to  be  optimal  are  not  satisfied, 
which  may  lead  to  excessive  false  alarms  and  missed  detections.  The  AOGF  and  OD 
are  based  on  noise  characteristics  that  better  describe  the  real  situation. 

We  shall  here  discuss  separately  the  implementation  of  two  algorithms;  adaptive  statisti¬ 
cally  optimal  group  filter  (AOGF),  here  also  refened  to  as  adaptive  beamformer  (AB),  and 
statistically  optimal  scalar  detector  (OD),  having  in  mind  that  the  two  will  be  combined  as 
a  two-step  detector  system. 

Time  consumption  for  the  detector  based  on  AOGF.  Broad  hand  AB  deployment. 

The  AOGF  procedure  consists  of  two  steps:  adaptation  and  filtering.  Adaptation  is  the 
most  time  consuming  and  is  not  realizable  in  real  time  (without  changing  to  a  recursive 
adaptive  beamformer). 

To  compute  the  covariance  matrix  function  and  the  vector  of  filter  coefficients  now  takes 
approximately  5  minutes  for  2  minutes  of  broadband  (40  Hz  sampling  rate)  25-channel 
data.  On  the  other  hand  AOGF  filtering  is  rather  fast  and  takes  approximately  the  same 
time  as  filtering  with  Butterworth  bandpass  filters. 

Since  the  procedure  is  broken  into  two  individual  steps,  the  adaptation  may  be  performed 
in  parallel  to  the  AOGF,  i.e.,  one  process  may  do  noise  adaptation  (generate  the  filter)  and 
report  this  into  a  data  base  that  may  be  used  by  the  filter  process.  So  the  AOGF  processing 
depends  only  on  how  often  the  parameters  of  the  group  filter  are  to  be  changed. 

In  Kushnir  et  al  (1990).  it  was  shown  that  high  noise  suppression  of  AOGF  without  read¬ 
aptation  remains  during  at  least  40  minutes.  On  Fig.  7.3.1  the  result  is  shown  of  AOGF 
noise  suppression  without  readaptation  for  50  hours.  Each  curve  deptcted  there  is  the  ratio 
of  noise  powers  P(average  channel  )/P(A(X;F)  and  P(beam)/P(AOGF).  By  beam,  we  mean 
the  conventional  array  beam  -  unfiltered. 
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In  this  noise  study  90  different  noise  time  intervals,  each  of  2  minute  length,  were  selected 
by  scanning  the  detection  list  and  requiring  that  the  noise  interval  started  at  least  6  minutes 
after  any  detection,  and  at  least  6  minutes  before  any  detection. 

Adaptation  was  performed  for  the  2  minutes  in  the  beginning  of  this  total  time  interval 
with  the  following  parameters:  frequency  band  0.2- 1.5  Hz,  resampling  factor  4,  azimuth 
334  degrees,  velocity  24  km/sec.  Subsequently,  AOGF  filtering  was  done  for  these  90 
intervals  of  noise. 

Note  that  the  resampling  factor  4  means  resampling  from  40  to  10  Hz,  i.e.,  the  adaptation 
is  here  done  for  the  frequency  range  0.2  -  1.5  Hz,  but  the  output  AOGF  and  conventional 
beam  represent  the  full  frec|uency  range  up  to  the  Nyquist  frequency  (5  Hz). 

From  Fig.  7.3.1  it  is  seen  that,  without  readaptation,  we  get  more  than  40  times  (16  dB) 
noise  suppression  compared  to  an  average  of  single  channels. 

Some  drops  in  the  ncuse  suppression  seen  on  the  figure  may  be  related  to  incoming  signals 
not  reported  in  the  detection  bulletin,  or  problems  with  finding  a  'clean'  noise  interval. 

The  similar  cuiacs  were  obtained  for  ARCESS  noise  with  noise  adaptation  to  the  fre¬ 
quency  bands  0  -  4.5  Hz  and  0.2  -  1.5  Hz  (Fig.  7.3.2  and  7.3.3). 

From  the  figures  it  is  seen  that  noise  suppression  by  AOGF  is  not  so  stable  as  for  NORESS 
but  nevertheless  high  noise  suppression  of  more  than  16  dB  is  consistently  achieved  for 
more  than  6  hours  after  adaptation. 

The  st.ibility  of  .-XOGF  noise  suppression  in  time  may  be  explained  by  the  fact  that  the 
.idaptation  is  done  in  the  lower  frequency  bands  w  here  noise  power  fluctuations  are  small¬ 
est  ( Fyen.  1 990). 

Let  us  consider  now  the  possibilities  of  constructing  a  small  aiTay  detector  system  consi.st- 
ing  of  several  adaptive  beams  steered  at  certain  azimuths  and  velocities.  Such  a  system 
w  ill  be  anakigous  to  the  on-line  NORESS  detector  system. 

The  main  question  here  is:  w  hat  is  the  minimum  number  of  AB  providing  adequate  detect¬ 
ability  while  still  allow  ing  real-time  or  routine  processing. 

To  ti7  to  answer  this  question,  an  experiment  was  performed  with  simulating  events  of  dif¬ 
ferent  SNR  -  signal  to  noise  ratio.  A  strong  regional  event  with  Pn  parameters:  azimuth 
193.  velocity  8.2.  frequency  band  2.5-4.5  Hz,  SNR  148,  was  selected.  This  event  was 
treated  as  a  signal  without  noise,  and  then  scaled  down  and  added  to  arbitrary  NORESS 
noise.  Figures  7.3.4  and  7.3.5  show  the  results.  The  5  curves  in  each  figure  correspond 
‘events'  of  5  different  SNR  as  function  of  slowness  of  the  AOGF  beam.  For  each  chosen 
azimuth  and  veltKity .adaptation  was  done  for  120  seconds  of  preceding  pure  noise,  fol¬ 
lowed  by  filtering  of  the  signal  segment  by  ACXJF  (i.e.,  forming  an  AOGF  beam).  The 
curves  show  SNR  for  the  AOGF  beam  (SNRCXJF)  and  the  corresponding  single  channel 
AO  maximum  SNR  for  each  different  scale  of  the  signal.  SNR  was  estimated  by  the  ratio 
of  STA  of  the  2  first  seconds  of  the  signal  and  average  STA  of  40  sec  of  pure  noise. 
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From  the  figures  we  see  that  we  get  highest  SNR  of  the  AOGF  at  the  correct  azimuth, 
whereas  SNR  is  nearly  constant  for  velocities  above  6  km/sec.  Looking  at  the  shape  of  the 
SNR  curve  for  azimuth,  we  see  for  azimuth  range  150  -  225  degrees,  the  SNR  varies  by  a 
factor  of  two. 

Allowing  a  predefined  beam  set  to  have  a  worst  case  missteering  loss  of  3  dB,  we  may 
infer  that  8  different  azimuths  and  two  velocities  may  be  enough  to  produce  high  detect¬ 
ability  for  all  azimuths  and  velocities.  That  is,  16  noise  adaptations  has  to  be  performed  at 
regular  intervals  separated  up  to  6  hours.  This  will  be  an  acceptable  load  for  most  systems. 


Frequency  dependence  of  noise  suppression  for  AOGF  and  standard  beams.  Low  fre¬ 
quency  AOGF  deployment. 

The  current  NORESS  detector  system  was  designed  mainly  for  the  detection  of  local  and 
regional  events,  having  its  best  performance  in  the  high  frequency  band  (>2  Hz).  But  as  it 
was  shown  by  Ringdal  (1990),  it  is  very  effective  also  for  the  detection  and  reporting  of 
teleseismic  signals. 

Most  of  the  teleseismic  events  detected  at  NORESS  are  from  the  Eurasian  continent  hav¬ 
ing  the  best  properties  of  propagation  at  high  frequencies  (1.5-4  Hz).  But  for  teleseismic 
signals  coming  from  other  directions,  energy  at  these  frequencies  attenuate  much  more 
strongly.  As  a  result,  much  of  the  teleseismic  information  is  masked  by  low  frequency 
noise  and  this  reduces  the  detectability  of  the  current  operational  system. 

We  shall  try  to  show  that  complementing  the  current  detector  system  with  SOA  may 
improve  the  detectability  also  for  teleseismic  distances.  Let  us  compare  noise  suppression 
by  AOGF  and  standard  beams  in  different  frequency  bands  reflected  in  Table  7.3. 1 .  By  a 
standard  beam  we  mean  one  of  the  conventional  beams  used  in  NORESS/ ARCESS  on¬ 
line  system.  Such  a  standard  beam  is  defined  by  array  sub-configuration,  velocity  and  azi¬ 
muth.  Sub-configurations  are  denoted  in  such  a  way  that,  e.g.,  AOAC,  means  the  center 
instrument  AO,  plus  the  A-ring  and  C-ring. 

The  table  comprises  STA  values  evaluated  on  pure  NORESS  noise  using  AOGF  and  stan¬ 
dard  beam  traces  filtered  in  the  most  important  frequency  bands,  that  is,  the  AOGF  as  the 
result  of  the  group  filtering  is  additionally  filtered  with  the  same  bandpass  filters  as  the 
suindard  beams.  The  AOGF  was  ba.sed  on  noise  adaptation  in  the  frequency  band  0-16  Hz. 
As  STA  is  an  estimate  of  the  power,  the  table  is  also  indicative  of  the  absolute  level  of 
noise  suppression. 

For  ACXjF.  the  25  vertical  NORESS  channels  (AOABCD  configuration)  were  used.  Noise 
adaptation  was  performed  using  a  100-second  noise  segment  preceding  the  40  seconds  of 
data  used  for  beams  and  AOGF  filtering.  For  these  teleseismic  AOGF  and  standard  beams, 
infinite  velocity  and  zero  azimuth  was  u.sed. 

From  the  table  it  is  seen  that; 

1 .  In  all  of  the  listed  frequency  bands,  the  resulting  noise  power  of  AOGF  is  lower  than 
the  standard  beams,  i.e.,  better  noise  suppression  is  obtained. 
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2.  The  best  standard  beam  configurations  were  AOD  for  0.3-1 .0  Hz  and  0.5-1.5  Hz,  and 
AOCD  for  1. 5-5.0  Hz  and  5-16  Hz. 

3.  Compared  to  the  AOD  beam,  broadband  AOGF  has  10  dB  better  suppression  in  the 
0.3-1 .0  Hz  and  6  dB  improvement  in  the  0.5-1. 5  Hz  frequency  bands. 

4.  In  the  highest  frequency  band  the  AOGF  and  AOCD  have  comparable  noise  suppres¬ 
sion  results. 

5.  As  the  higheut  frequency  of  the  lowpass  prefiltering  in  the  noise  adaptation  (lowpass 
16  Hz  in  the  first  part  of  the  table)  decreases,  the  STA  power  of  AOGF  also 
decreases,  that  is.  noise  suppression  improves. 

For  a  lowpass  (Ip)  1.5  Hz  filter,  the  AOGF  shows  approximately  6  dB  smaller  STA  value 
for  the  0.3- 1 .0  Hz  and  0.5- 1 .5  Hz  bands,  as  opposed  to  adaptation  with  Ip  16  Hz.  In  the 
band  0.2-0.6  Hz,  comprising  a  considerable  part  of  the  noise  power,  the  AOGF  is  sup¬ 
pressed  1 9  dB  compared  to  AOD. 

Hence,  the  best  teleseismic  perfonnance  of  AOGF  is  achieved  by  prefiltering  the  data  with 
a  low  pass  filter  1.5  Hz 


Comparison  of  optimal  detector  and  STA/LTA  detectability 

We  will  here  present  some  preliminai7  results  concerning  application  of  the  one-dimen¬ 
sional  statistically  optimal  detector,  i.e.,  the  so-called  ESTDET  algorithm  (Kushnir  et  al, 
1991 ).  This  algorithm  is  based  on  optimal  accounting  of  noise  features.  Compared  to  the 
ordinary  STA/LTA  algorithm,  this  algorithm  not  only  takes  into  account  amplitude  level, 
but  also  variations  of  the  spectrum  content. 

The  pnKedure  has  two  steps:  The  first  is  adaptation  on  a  time  segment  of  pure  noise.  As  a 
result  we  get  an  estimate  of  whitening  filter  coefficients.  The  second  step  is  filtering  using 
the  coefficients  computed  at  the  previous  step  and  computation  of  Chi-square  statistic  in  a 
moving  window. 

To  examine  the  performance  of  ESTDET  as  a  function  of  SNR,  we  have  chosen  a  NOR- 
ESS  recording  of  a  regional  event  in  Sweden,  i.e.,  the  same  event  used  for  Fig.  7.3.4. 

The  event  was  detected  with  relatively  high  SNR  by  standard  beam  N055  (configuration 
AOBCD,  frequency  band  2.5-4.5  Hz).  Again  we  used  the  procedure  of  superposing  this 
signal  onto  pure  NORESS  noise  with  different  .scaling. 

Fig.  7.3.6  shows  the  display  of  the  bandpass  filtered  beam  N055,  the  STA  for  this  beam, 
and  ESTDET  (XlSQTl  Chi-square  statistics  from  ESTDET  process)  for  this  original  sig¬ 
nal.  The  lower  part  of  the  figure  shows  the  same  three  processes  when  the  signal  is  scaled 
down  by  a  factor  80  in  amplitude  and  added  to  a  noise  trace. 

The  advantage  of  ESTDET  over  STA  for  the  latter  case  is  obvious. 
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The  next  experiment  was  to  obtain  ESTDET  statistics  with  varying  noise  scale  factors  for 
the  simulated  data. 

The  ESTDET  statistic  was  computed  for  the  AOGF  beam  and  for  the  standard  beam.  Sig- 
nal-to-noise  ratio  was  estimated  by  the  ratio  of  ESTDET  for  the  3  first  seconds  of  P-wave 
and  40  seconds  of  preceding  pure  noise,  i.e.  Xi-square  values  are  averaged  for  the  noise 
segment  to  get  a  signal-to-noise  ratio  comparable  to  STA/LTA.  For  the  standard  beam,  the 
STA  was  used. 

Fig.  7.3.7  shows  the  results.  The  upper  curve  is  the  ratio  of  ESTDET  statistic  maxima  on 
signal  and  noise.  The  curve  is  computed  for  unfiltered  AOGF  data.  The  bottom  curve  is 
the  corresponding  ratio  for  the  STA  algorithm  (SNR)  on  standard  beam  N055.  This  ratio 
estimates  STA/LTA.  The  middle  curve  is  similar  to  the  upper  one,  but  additionally  band¬ 
pass  filtered  in  the  frequency  band  2.5-4.5  Hz. 

Comparing  upper  and  bottom  curves,  it  is  seen  that  the  ESTDET  SNR  is  better  than  that  of 
the  STA/LTA  algorithm,  especially  for  the  unfiltered  AOGF  beam.  For  example,  for  noise 
scale  equal  to  100,  SNR  for  STA  is  2  and  for  ESTDET  is  50-60. 

To  examine  ESTDET  false  alarms  we  have  chosen  from  two  days  bulletin  all  noise  inter¬ 
vals  of  at  least  2  minutes  length  and  selected  2  minutes  segments  from  each.  The  total 
amount  of  noise  is  about  3  hours. 

Fig.  7.3.8  shows  the  number  of  ESTDET  values  exceeding  a  given  threshold,  as  a  function 
of  threshold.  For  low  thresholds  it  is  very  lai^e,  but  it  decreases  with  increasing  of  the 
threshold  value.  Only  one  exceedance  is  detected  for  a  threshold  larger  than  5()0.  Let  us 
compare  it  with  the  maximum  values  of  ESTDET  on  the  signal+noise  (Fig.  7.3.9)  com¬ 
puted  for  different  noise  scales.  It  is  seen  that  value  5(K)  corresponds  to  a  noise  scale  of  85 
or  STA/LTA  value  around  2  which  is  lower  than  the  conventionally  used  threshold. 

So,  we  may  suppose  that  if  the  threshold  for  ESTDET  is  equal  to,  say,  550,  it  is  possible  to 
detect  a  weak  regional  P-wave,  STA/LTA=2,  with  a  low  probability  of  false  alarms  during 
several  hours  of  processing. 

From  these  preliminary  results  we  may  draw  the  following  conclusions; 

1 .  AOGF  is  computed  after  adaptation  to  the  noise.  This  process  is  time  consuming,  but 
it  is  shown  that  high  noise  suppression  may  be  obtained  using  only  one  noise  adapta¬ 
tion  for  several  hours  of  data  processing. 

2.  It  has  been  shown  that  for  a  NORESS-type  array  adaptation  to  the  noise  for  about  8 
different  azimuths  and  two  velocities,  is  sufficient  to  obtain  adequate  beam  coverage 
for  regional  P  phases. 

3.  The  main  advantage  of  AOGF  is  to  be  expected  for  low  frequency  teleseismic  sig¬ 
nals,  where  we  have  both  good  detectability  and  a  broadband  undisturbed  estimate 
of  the  signal. 
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4.  The  ESTDET  detector  is  much  more  sensitive  to  the  detection  of  seismic  phases 
compared  to  conventional  STA/LTA  due  to  optimal  accounting  to  noise  spectrum 
variations. 

5.  Due  to  higher  sensitivity  the  ESTDET  process  may  be  used  for  the  detection  and 
timing  of  much  weaker  seismic  phases  than  detected  now. 

V.  Pinsky,  MITPAN  Institute,  Moscow 
S.  Tsvang,  MITPAN  Institute,  Moscow 
J.  Fyen 
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Adaptation  Ip  16  Hz 


Filter  band 

AOGF 

AOABCD 

AOBC 

AOD 

AOCD 

AOAB 

AO 

0.3  - 1.0  Hz 

0.035 

0.150 

0.188 

0.’05 

0.131 

0.206 

0.207 

0.5-  1.5  Hz 

0.044 

0.142 

0.192 

0.085 

0.114 

0.219 

0.226 

1.5  -5.0  Hz 

0.018 

0.038 

0.062 

0.032 

0.022 

0.093 

0.113 

5.0  -  16.0  Hz 

0.005 

0.0056 

0.0076 

0.065 

0  0071 

0.0094 

0.027 

Unfiltered 

0.05 

0.174 

0.225 

0.157 

0.145 

0.256 

0.271 

Adaptation  Ip  4.5  factor  4 

0.3  -  1 .0  Hz  0.023 

0.5 -1.5  Hz  0.029 

1.5 -5.0  Hz  0.016 

Unfiltered  0.034 

Adaptation  Ip  1.5  factor  4 

0.2 -0.6  Hz  0.011  0.090  0.139 

0.3-1.  Hz  0.019 

0.5 -1.5  Hz  0.023 

Untillercd  0.025 

Table  7J.1.  Comparison  of  AOGF  noise  suppression  in  different  frequency  bands.Ca.se 
study  of  NORESS  noise  segment  start  1991-140:00.05.00.004.  The  table  shows  STA 
power  values  (normalized  with  factor  KXX))  for  AOGF  and  different  standard  beams. 
■AOGF  adaptation  using  velocity  99999,  azimuth  0.0.  lowpass  16  Hz  prefilter  and  100  sec¬ 
onds  noise.  Data  are  not  resampled  (resample  factor  1).  STA  power  values  on  the  follow¬ 
ing  40  seconds.  Before  STA,  traces  arc  butterworth  bandpass  filtered  in  the  given  filter 
bands.  AOBC.  etc.,  show  which  sub-configuration  that  is  used.  Additional  values  are  com¬ 
puted  for  AOGF  using  lowpass  4.5  and  1 .5  Hz  and  resampling  factor  4.  (Resample  from 
40  to  in  Hz). 


90 


GAINOGF 


GAINOGF 


NORSAR  Sci.  Rep.  1-91/92 


November  1991 


OGF  GAIN  VERSUS  TIME  WITHOUT  READAPTATION  (0--4.5Hz) 


start  1991-174:12.20  TIME(hours) 


Fig.  7J.2.  AOGF  gain  versus  time  without  readaptation  on  pure  ARCESS  noise  selected 
from  bulletin.  OGFGAIN  is  the  ratio  of  noise  powers  on  2-minute  noise  intervals.  Start 
time  1991-174:12.20  Adaptation  interval  1991-176:13.20  -  13.22.  Data  frequency  band  0- 
4.5  Hz.  Adaptation  velocity  10  km/sec.  azimuth  193  degrees. 
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OGF  GAIN  VERSUS  TIME  WITHOUT  READAPTATION  (0.2-1. 5  Hz) 


15  20  25  30 

start  1991-174:12.20  TIME(hours) 


Fig.  733.  AOGF  gain  versus  time  without  readaptation  on  pure  ARCESS  noise  selected 
from  bulletin.  OGFGAIN  is  the  ratio  of  noise  powers  on  2-minute  noise  intervals.  Start 
time  1991-174:12.20  Adaptation  interval  1991-176:13.20  - 13.22.  Data  frequency  band 
0.2- 1.5  Hz.  Adaptation  velocity  10  km/sec,  azimuth  193  degrees. 
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OGF  AZIMUTH  MIS-STEERING 


AZIMUTH  (true  azimuth  193  deg.) 


Fig.  7  J.4.  AOGF  azimuth  missteering  on  regional  event  NORESS  P-wave  with  different 
scaling  of  added  noise.  Signal  parameters;  start  1991-086:05.18.46,  velocity  8.2  km/sec 
azimuth  193  degr^,  freq.  3.35  Hz.  Noise  start  1991-140:00.05.  OAGF  adaptation  on  100 
seconds  of  pure  noise  with  velocity  8.2  km/sec,  frequency  band:  0-4.5  Hz,  and  different 
azimuths.  Filtering  on  50  .seconds  of  signal  +  noise.  SNR(X}F=  (average  STA  on  2  sec.  of 
P-wavc)/{average  STA  on  40  sec.  of  noise) 
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OGF  VELOCITY  MIS-STEERING 


VELOCITY  (true  vel=8.2  km/sec) 


Fig.  7  J.5.  AOGF  velocity  missteering  on  regional  event  NORESS  P-wave  with  different 
scaling  of  added  noise.  Signal  parameters;  start  1991-086:05.18.46,  velocity  8.2  km/sec 
azimuth  193  degrees,  freq.  3.35  Hz.  Noise  .stan  1991-140:00.05.  OAGF  adaptation  on  100 
seconds  of  pure  noise  with  azimuth  193  degrees,  frequency  band:  0-4.5  Hz  and  different 
velocities.  Filtering  on  50  seconds  of  signal  +  noise.  SNRC)GF=  (average  STA  on  2  sec.  of 
P-wave )/(average  STA  on  40  sec.  of  noise) 


95 


NORS AR  Sci.  Rep.  I  -9 1/92 


November  1991 


'99'  - C  S  6 
sc  c  i  e  —  SC 


N.-'  •  a/9  •  ’  1  jg.a3 


Fig.  73.6.  STA  for  ESTDET  and  standard  beam  n055  for  regional  event.  Signal  parame¬ 
ters:  start  1991-086:05.18.46.  velocity  8.2  km/sec  azimuth  193  degrees,  freq.  3.35  Hz. 
Noise  start  1991-140:00.05.  Traces:  1.  standard  bean  (SB),  ff.band  2.5-4.5  Hz;  2.  STA  on 
SB;  3.  ESTDET  on  SB;  4.  SB  +  noise,scale=80;  5.,6.  STA  and  ESTDET  on  SB  +  noise. 
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ESTDET(w=1,o=3,n=10):  SIGNAL(dlst=258,v=8.2,az=l  92) 


SNRestdet  =  MAXsIgnal/MAXnolse 

*  ESTDET  after  unfiltered  beam 
A  ESTDET  after  filtered  beam 
o  STA/LTA  after  standard  beam 


Fig.  73.7.  Comparison  of  ESTDET  and  STA/LTA  on  the  regional  P-phase  +  noise  for  the 
different  noise  scales.  ESTDET  parameters;  (w=l,o=3,n=l0). 
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ESTDET(w=  1  .o=3,n  =  1 0)  SIGNAL(clist=258,v=8.2,az=  1  92) 


Fig.  7.3.9.  Maximum  values  of  estdet  on  regional  P-phase  +  noise  versus  noise  scale 
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7.4  Extensions  of  the  Northern  Europe  Regional  Array  Network  —  a 
new  three-component  station  at  Apatity,  USSR,  and  a  planned 
array  at  Spitsbergen 

The  Northern  Europe  Regional  AiTay  Network  currently  comprises  the  high-frequency 
arrays  NORESS  and  ARCESS  in  Norway,  FINESA  in  Finland  and  GERESS  n  Germany, 
as  well  as  the  two  three-component  stations  KSP  and  SFP  in  Poland.  The  locations  of 
these  stations  are  shown  in  Fig.  7.4.1  Data  from  all  of  these  installations  are  transmitted 
continuously  and  in  real  time  to  the  NORSAR  Data  Processing  Center  at  Kjeller,  Norway, 
where  the  data  are  processed  jointly  by  the  Intelligent  Monitoring  System  (IMS).  Current 
prospects  for  extending  this  network  and,  if  possible,  integrating  new  data  .sources  into 
IMS,  are  described  in  this  contribution. 

The  establishment  of  the  FINESA  luray  in  Finland  was  made  possible  through  joint  efforts 
between  NORSAR  and  the  University  of  Helsinki,  based  on  a  bilateral  cooperative  agree¬ 
ment.  A  similar  arrangement  has  recently  been  made  with  the  Kola  Science  Center  of  the 
USSR  Academy  of  Sciences,  and  as  first  result  of  this  cooperation,  a  high-quality  three- 
component  station  was  installed  by  NORSAR  personnel  in  Apatity  (see  Fig.  7.4.1 )  in  June 
1991.  The  following  paragraph  gives  a  description  of  the  system  installed  in  Apatity  and 
contains  material  related  to  analysis  of  data  from  this  new  station. 

According  to  current  plans,  a  9-element  array  with  a  1 -kilometer  aperture  will  be  installed 
on  the  island  of  Spitsbergen  in  the  Arctic  (see  Fig.  7.4.1 )  in  the  summer  of  1992.  A  short 
description  of  these  plans  are  given  in  this  contribution. 

A  new  3 -component  station  at  Apatity.  USSR 

In  mid-June  1991,  a  high-quality  three-component  .station  was  installed  in  Apatity,  on  the 
Kola  Peninsula  of  the  USSR.  The  station  was  installed  and  is  operated  under  an  agreement 
on  scientific  cooperation  within  seismology  between  the  Kola  Science  Center  of  the  USSR 
Academy  of  Sciences  and  NORSAR. 

The  three-component  station  installed  in  Apatity  comprises  S-13  .seismometers,  a  Nano¬ 
metrics  RD3  digitizing/synchronization/multiplexing  unit,  an  Omega  clock  and  a  PC- 
based  system  for  data  acquisition,  analysis  and  archiving.  The  amplitude  response  of  the 
system  is  shown  in  Fig.  7.4.2.  The  PC  system  has  a  25  MHz  386  processor,  4  Mbyte 
RAM.  a  300  Mbyte  ESDI  disk,  and  a  cartridge  tape  drive.  The  system  was  installed  in  the 
basement  of  the  building  of  the  Seismological  Laboratory  of  the  Kola  Science  Center,  in 
the  town  of  Apatity.  The  S-13  seismometers  were  placed  on  the  same  pad  as  the  seismom¬ 
eters  of  the  analog  APA  station  (established  in  1956).  Data  are  recorded  continuously  on 
magnetic  disk  in  Apatity  and  copied  to  cartridge  tapes,  which  are  sent  to  NORSAR.  An 
archive  of  continuous  data  from  the  new  station  in  Apatity  is  thus  maintained  both  at  Apa¬ 
tity  and  at  NORSAR. 
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Analysis  of  data  recorded  at  the  new  Apatitv  station  --  noise  levels 

Fig.  7.4.3a  shows  six  noise  spectra  for  the  new  Apatity  station.  All  noise  samples  are  taken 
during  daytime  (2  pm  to  5  pm  local  time),  with  the  exception  of  the  sample  corresponding 
to  the  dashed  curve,  which  is  taken  around  midnight,  local  time.  The  trend  of  the  curves 
strongly  indicates  cultural  origin  for  the  noise  above  2  Hz.  Also,  one  must  keep  in  mind 
that  the  station  is  located  within  the  town  of  Apatity,  which  has  approximately  100,000 
inhabitants.  The  very  heavy  mining  activities  in  the  nearby  (20  km)  Khibiny  Massif  is  sup¬ 
posedly  also  contributing  to  the  noise  level.  A  distinct  and  consistent  peak  at  4  Hz  indi¬ 
cates  a  localized  source.  This  4  Hz  wave  can  be  clearly  seen  in  the  noise  preceding  the 
onset  in  the  recording  of  the  shot  in  Mine  II  in  Fig.  7.4.5.  All  spectra  except  one  also 
exhibit  a  pronounced  peak  at  2.5  Hz.  The  origin  of  these  peaks  in  the  noise  spectra  are  as 
yet  unknown.  There  is  an  indication  of  a  slightly  lower  noise  level  above  7-8  Hz  at  night¬ 
time. 

Compared  to  the  NORESS  noise  model  (superimposed  on  Fig.  7.4.3a),  the  noise  level  at 
the  new  Apatity  station  .stands  out  as  being  very  high  above  2  Hz.  For  lower  frequencies, 
the  noise  levels  at  NORESS  and  Apatity  are  similar.  This  is  also  to  be  expected  on  the 
basis  of  usual  assumptions  about  noise  generating  mechanisms  for  frequencies  below  2 
Hz. 

On  12  June,  prior  to  installing  the  station  at  the  Seismological  Laboratory  in  Apatity,  the 
entire  system  was  deployed  for  an  hour  at  a  location  approximately  15  km  to  the  west  of 
Apatity  (denoted  PROF  in  Fig.  7.4.4).  This  is  a  candidate  location  for  a  possible  future 
deployment  of  a  small  array  in  Apatity.  While  the  noise  level  is  still  high  relative  to  the 
NORESS  noise  model,  the  level  is  clearly  lower  than  in  the  town  of  Apatity,  as  can  be 
.seen  by  comparing  Fig.s.  7.4.3a  and  7.4.3b. 

Analysis  of  data  recorded  at  the  new  Apatitv  station  -  P-wave  signals 

The  map  of  Fig.  7.4.4  shows  the  location  of  the  new  Apatity  station  and  the  nearby  mining 
area  of  the  Khibiny  Massif.  The  lateral  extent  of  the  mining  area  is  approximately  15  km, 
and  there  is  information  available  that  allows  the  assignment  of  each  individual  explosion 
since  early  1991  to  one  of  the  five  subareas  indicated  in  Fig.  7.4.4.  Fig.  7.4.5  shows  three- 
component  records  for  the  new  Apatity  station  for  five  such  explosions,  one  from  each  of 
the  five  mines.  Detailed  information  on  these  shots  is  given  in  Table  7.4.1.  From  Fig. 
7.4.5.  it  is  clearly  seen  that  there  is  a  lai^e  variation  among  the  shots  with  respect  to  ampli¬ 
tude  and  frequency  content  of  P,  S  and  Rg  waves,  and  this  variation  cannot  be  explained 
by  differences  in  epicentral  distance  and  shot  size  alone.  Hence,  other  factors  like  propa¬ 
gation  path  effects  and  specific  conditions  pertaining  to  the  individual  shots  are  important. 
For  example,  information  on  shot  geometry  and  ripple  firing  delay  times  is  not  currently 
available. 

As  a  check  on  the  performance  and  integrity  of  the  new  three-component  Apatity  station, 
we  used  the  data  for  the  five  shots  in  Fig.  7.4.5  to  derive  the  apparent  arrival  azimuth  for 
the  P  wave  onset  using  the  broad  band  slowness  estimation  technique  of  Kviema  and 
Doombos  (1986).  The  results  are  given  in  Table  7.4.1 .  along  with  the  true  azimuths,  which 
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are  given  in  terms  of  intervals,  due  to  the  spatial  extent  of  each  of  the  five  mines.  The  cal¬ 
culations  were  made  for  four  different  frequency  bands,  as  shown  in  the  table.  The  length 
of  the  time  windows  used  were  1.0  s  for  the  8-16  Hz  band,  1.5  s  for  the  5-10  Hz  and  3-5 
Hz  bands,  and  2.0  s  for  the  1-3  Hz  band.  The  start  time  for  all  windows  was  0.1  s  prior  to 
the  P  arrival  time.  The  table  shows  reasonable  figures  with  an  average  deviation  of  5.5 
degrees  from  the  midpoint  of  the  interval  for  the  true  azimuth,  when  disregarding  the  low 
SNR  event  from  mine  II.  We  take  this  result  to  indicate  that  the  new  Apatity  station  per¬ 
forms  well  from  a  technical  point  of  view. 

As  a  further  illustration  of  the  new  data  now  available  from  Apatity,  recordings  at  Apatity, 
ARCESS  and  FINESA  of  a  probable  earthquake  approximately  200  km  due  south  of  Apa¬ 
tity  are  shown  in  Fig.  7.4.6.  The  bulletin  of  the  University  of  Helsinki,  Finland,  gives  an 
epicenter  of  65.71®N,  33.24®E,  based  on  altogether  41  readings  of  P  and  S  phases.  The 
magnitude  is  Ml  3.7  and  the  depth  is  estimated  at  8.8  km  (standard  deviation  4.5  km)  in 
the  Helsinki  bulletin.  The  arrival  azimuths  of  the  first  P  wave  at  the  new  Apatity  station 
were  estimated  as  was  done  for  the  five  mining  shots,  and  results  are  also  given  in  Table 
7.4.1 .  Again,  the  results  are  reasonable,  with  an  average  deviation  from  the  true  azimuth  of 
about  6  degrees.  A  feature  of  special  interest  in  Fig.  7.4.6  is  the  Rg  wave  seen  in  the  Apa¬ 
tity  and  FINESA  data.  Especially  surprising  is  the  presence  of  an  Rg  phase  at  FINESA,  at 
a  distance  of  593  km.  The  strong  Rg  waves  also  point  toward  a  shallow-focus  event,  and 
maybe  shallower  than  8.8  km,  as  given  in  the  Helsinki  bulletin. 

We  think  that  the  new  data  now  available  from  Apatity  will  prove  to  be  very  valuable.  It  is 
of  special  interest  that  we  are  now  able  to  obtain  near-field  recordings  of  the  mining  shots 
in  the  Khibiny  Massif.  All  of  these  shots  are  well  recorded  at  ARCESS  (distance  480  km), 
and  the  majority  of  the  larger  ones  also  at  FINESA  (distance  780  km)  and  NORESS  (dis¬ 
tance  1315  km).  Plans  including  a  computer-to-computer  link  between  NORSAR  and  the 
Kola  Science  Center  and  possible  deployment  of  additional  stations  are  row  being  dis- 
cus.sed. 


Plans  for  a  small-aperture  array  at  Spitsbergen 

Companies  taking  part  in  oil  exploration  and  production  on  the  Norwegian  Continental 
Shelf  have  earlier  sponsored  the  establishment  of  a  network  of  seismic  stations  onshore 
northern  Norway.  As  a  further  development  of  this  network,  NORSAR  will  establish, 
under  the  sponsorship  of  the  same  oil  companies,  a  small-aperture  array  near  Longyear- 
byen  on  the  island  of  Spitsbergen  during  the  summer  of  1992. 

A  site  survey  was  conducted  in  August  1991 .  A  suitable  location  was  found  approximately 
20  km  from  Longyearbyen,  which  is  the  largest  Norwegian  settlement  and  also  adminis¬ 
trative  center  on  the  island.  The  plan  is  to  install  a  9-element  1 -kilometer  aperture  array  of 
vertical  sensors.  The  array  will  have  two  concentric  rings  with  3  and  5  elements,  respec¬ 
tively.  plus  one  sensor  at  the  center.  The  rocks  at  the  site  are  of  Cretaceous  age,  covered  by 
thin  moraine  of  variable  depth.  The  sensors  will  be  placed  in  boreholes,  drilled  either  to 
the  Cretaceous  rock  or  to  a  depth  where  stable  permafrost  conditions  are  obtained  (depths 
of  the  order  of  6  meters  are  considered  to  be  sufficient  in  this  regard).  All  necessary  per¬ 
missions  to  deploy  the  array  have  been  granted  by  the  appropriate  authorities. 
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The  site  survey  included  noise  measurements  at  the  .selected  site.  A  noise  spectrum  is 
shown  in  Fig.  7.4.7.  Even  if  the  noise  level  is  seen  to  be  higher  than  the  average  noise 
level  at  NORESS,  we  consider  the  noi.se  level  at  the  Spitsbeigen  site  to  be  quite  satisfac¬ 
tory  for  an  array  in.stallation.  For  example,  there  are  no  specific  peaks  in  the  spectrum. 
There  is  mining  activity  in  the  area,  with  mines  located  at  distances  between  10  and  20  km 
from  the  site. 

The  system  to  be  established  at  Spitsbergen  will  have  an  on-site  recording  capability  only. 
We  do,  however,  consider  that  the  data  from  this  new  array  would  be  of  great  value  to  the 
IMS.  This  requires  continuous  transmission  via  satellite  of  data  from  Spitsbergen  to  Kjel- 
ler.  The  annual  lease  of  such  a  link  is  estimated  at  approximately  USD  30,000. 

S.  Mykkeltveit  A.  Dahle  J.  Fyen  T.  Kvaerna 

P.W.  Larsen  R.  Paulsen  F.  Ringdal 

E.O.  Kremenetskaya,  Kola  Science  Center 
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Epicenter 

Shot  size 

P  wave 

True 

Estimated  P-wave  arrival  azimuth  (deg) 

(tons) 

arrival  time 

azimuth 

(deg) 

1-3  Hz 

3-5  Hz 

5-10  Hz 

8-16  Hz 

Mine  I 

120 

216/12.49.55.8 

36-46 

47.1 

38.5 

44.6 

37.0 

Mine  11 

35 

259/12.19.27.1 

46-53 

(55.3) 

(99.5) 

(78.9) 

(34.0) 

Mine  III 

94 

187/19.24.01.9 

55-60 

58.5 

57.6 

68.2 

71.7 

Mine  IV 

430 

207/12.31.00.3 

60-66 

69.2 

61.6 

74.2 

63.7 

Mine  V 

392 

193/10.21.22.1 

68-74 

72.4 

64.5 

60.1 

64.8 

65.71“N.33.24“E 

- 

236/10.57.01.2 

180.2 

174.3 

182.8 

185.5 

189.4 

Table  7.4.1.  Apparent  anival  azimuths  at  vhe  new  3-component  station  at  Apatity,  as 
determined  from  broad  band  slowness  estimation  using  the  method  of  Kvaema  and  Doom- 
bos  (1986).  The  epicentral  information  on  the  earthquake  at  the  bottom  of  the  table  is  from 
the  Helsinki  bulletin.  The  SNR  for  the  event  from  Mine  II  is  very  low  in  all  frequency 
bands  and  the  estimated  arrival  azimuths  are  given  in  parenthesis. 
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Fig.  7.4.1.  The  map  shows  the  existing  stations  of  the  Northern  Europe  Regional  Airay 
Network  as  well  as  the  location  of  the  two  new  stations  at  Apatity  and  Spitsbei:gen,  dealt 
with  in  this  contribution. 
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Fig.  7,4J.  a)  (top);  Six  noise  spectra  for  the  new  Apatity  station,  b)  (bottom);  Noise  spec¬ 
trum  from  the  location  outside  Apatity,  denoted  PROF  in  Fig.  7.4.4.  The  .solid  line  in  each 
of  the  figures  denotes  the  NORESS  noi.se  model  of  Fyen  (1990).  All  spectra  were  esti¬ 
mated  from  partly  overlapping  5-second  windows,  using  the  indirect  covariance  method 
(see  Fyen,  1990).  The  total  length  of  each  noise  sample  was  approximately  1  minute. 
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Fig.  7.4.4.  Map  of  the  Apatity  area,  .showing  the  location  of  the  new  Apatity  station  and 
also  the  candidate  site  PROF  for  a  future  small  array.  Also  shown  is  the  mining  area  of  the 
Khibiny  Massif,  and  the  division  into  five  (denoted  I  -  V)  separate  subareas  where  explo¬ 
sions  are  currently  taking  pla;e.  Locations  for  25  mining  induced  earthquakes  (magnitude 
range  2-5)  are  also  shown  (see  also  section  7.7,  this  volume.). 
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Mine  IV  Mine  V 


Fig.  7.4.5.  Apatity  three-component  data  for  five  mining  shots,  one  from  each  of  the  mines 
I  -  V  in  Fig.  7.4.4. 
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F^.  7.4.6.  Apatity,  ARCESS  and  FINESA  single  vertical  component  recordings  for  a  pre¬ 
sumed  earthquake  200  km  south  of  Apatity.  For  each  station,  the  plot  shov^rs  unfiltered 
data,  and  data  filtered  in  the  band  0.2-0.8  Hz,  Note  the  presence  of  an  Rg  phase  at  FINESA 
at  a  distance  of  nearly  600  km.  The  time  axis  refers  to  the  Apatity  station  on  top;  ARCESS 
and  FINESA  data  are  aligned  with  respect  to  the  P  onset. 
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Fig.  7.4.7.  Noise  spectrum  taken  at  the  site  for  the  new  small-aperture  array  to  be  installed 
on  Spitsbergen  in  1992.  The  solid  line  denotes  the  NORESS  noise  model  of  Fyen  (1990). 
The  spectrum  was  estimated  in  the  same  way  as  the  spectra  of  Fig.  7.4.3. 
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7.5  Initial  testing  of  mixed  event  separation  using  a  statistically 
optimal  adaptive  algorithm. 


Introduction 

The  networic  of  the  three  Fennoscandian  regional  arrays,  ARCESS,  FINESA  and  NOR- 
ESS,  has  shown  an  excellent  capability  to  monitor  the  seismic  activity  at  the  Soviet  test 
site  at  Novaya  Zemlya.  In  a  previous  report  (Kvaema  and  Ringdal,  1990),  we  demon¬ 
strated  the  application  of  the  threshold  monitoring  (TM)  technique  to  obtain  continuous 
estimates  of  the  largest  size  of  events  that  might  go  undetected  by  the  monitoring  network. 
This  study  showed  that  during  normal  noise  conditions,  it  was  possible  to  monitor  the 
Novaya  Zemlya  test  site  in  this  way  down  to  mj,  2.5  more  than  99%  of  the  time. 

However,  during  intervals  of  increased  seismic  amplitude  levels  at  the  arrays,  e.g.,  caused 
by  the  arrivals  of  signals  from  strong  earthquakes,  the  monitoring  capability  may  deterio¬ 
rate  significantly,  as  shown  in  Fig.  7.5.1.  For  the  day  analyzed  in  Fig.  7.5.1,  a  maximum 
network  threshold  of  3.8  occuned  because  of  airivals  from  a  strong  earthquake  at  Hindu 
Kush  (mj,  =  6.0).  This  implies  that  an  event  of  magnitude  3.8  in  theory  could  have 
occurred  at  Novaya  Zemlya  at  that  time  without  exceeding  the  network  threshold. 

To  improve  the  monitoring  capability  during  such  intervals  of  increased  magnitude  thresh¬ 
olds.  we  have  started  to  evaluate  the  potential  of  a  statistically  optimal  processing  scheme 
(Kushnir  and  Lapshin,  1984)  for  suppressing  the  effect  of  interfering  coherent  arrivals. 
The  method,  denoted  adaptive  optimal  group  filtering  (AOGF),  has  earlier  been  applied 
for  extraction  of  signals  from  background  noise  (Kushnir  et  al.,  1990).  Compared  to  con¬ 
ventional  beamforming  at  NORESS  and  ARCESS,  AOGF  has  demonstrated  a  capability 
to  provide  significant  SNR  gains  over  a  wide  frequency  band  (0.2-5.0  Hz).  The  largest 
gain  (12-18  dB  more  than  conventional  beamforming)  is  obtained  at  low  frequencies 
where  the  method  utilizes  the  very  coherent  nature  of  the  noise.  But  even  for  the  frequency 
band  from  2  to  4  Hz,  AOGF  provides  a  typical  SNR  gain  of  the  order  of  6  to  8  dB  com¬ 
pared  to  the  conventional  beam. 

Theoretical  background  for  adaptation  to  interfering  events 

In  the  study  of  Kushnir  et  al.  (1990),  a  pure  noise  interval  preceding  the  signal  was 
selected  for  adaptation.  Exploiting  the  stationarity  of  the  background  noise  field,  the  inter¬ 
val  containing  the  signal  was  subsequently  processed  by  the  AOGF  method  using  the 
parameters  obtained  during  the  adaptation. 

This  procedure  is.  however,  not  relevant  when  we  want  to  extract  a  hypothetical  signal 
hidden  somewhere  in  the  coda  of  another  signal.  In  this  case,  we  do  not  know  if  or  when 
the  hypothetical  event  occurred,  and  do  only  have  information  on  the  slowness  and  azi¬ 
muth  of  the  hypothetical  signal.  For  AOGF  to  be  applicable  to  our  new  purpose,  the  adap¬ 
tation  needs  to  be  successful  even  if  the  adaptation  interval  contains  the  wave  we  want  to 
extract. 
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We  will  in  the  following  show  theoretically  that  the  AOGF  method  has  the  desired  prop¬ 
erty: 

The  AOGF  method  is  based  on  an  optimal  Wiener  group  filter  (OGF)  with  frequency 
response: 


gif)  =  ig\if),-;gnif)) 


(1) 


where  m  is  the  number  of  array  sensors. 

The  OGF  has  the  characteristics  of  effectively  suppressing  the  coherent  part  of  the  noise 
field  while  keeping  the  signature  of  the  wanted  plane  wave  undistorted.  The  filter  ^  ^ 
transforms  the  input  array  data  =  (jc,,,  (where  T  denotes  the  transpose)  into  a 

scalar  trace  y, .  v,  =  |,*  .t,  where  *  denotes  convolution,  and  y(f)  =  g(J)  xij),  where  y  (J) , 
g  (f)  and  x(/)  are  the  Fourier  transforms  of  y,,  and  x,. 

It  is  well  known  that  |  (f)  has  the  form: 


gif)  =  -T - 

if)  ■r„^if)  -hif) 


(2) 


where  F,  (/)  is  the  power  spectral  density  matrix  of  the  noise  h,  =  (/i,,, 

h(J)  =  . (/) )  ^.  where  h^lf)  =  e  and  is  the  time  delay  at  the  *'*  sensor 

for  the  plane  wave  to  be  extracted  from  the  background  noise  field,  t  denotes  the  Henni- 
tian  conjugate. 

If  we  are  adapting  to  a  pure  noise  interval,  the  matrix  (/)  is  estimated  by  a  special  pro¬ 
cedure.  From  a  computational  and  statistical  point  of  view,a  good  way  to  estimate  the 
large  size  matrix  (/)  is  to  use  autoregressive  modelling  of  the  noise  field  n,.  The  new 
filter 


!(/)  = 


^'^if)-K'if) 
if) if) -h  if) 


(3) 


is  referred  to  as  the  adaptive  optimal  group  filter  (ACXjF)  where  ^„'  (/)  is  the  autoregres¬ 
sive  estimate  of  F^'  (/3 . 


If  the  adaptation  interval  contains  a  mixture  of  the  signal  and  pure  noise,  i.e.,  x,  =  s,  +  h„ 
we  cannot  estimate  F^'  O')  but  only  the  spectral  density  matrix  F^'  (J)  of  the  mixture.  Let 
us  consider  the  group  filter 
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d{f)  = 


(/)  •  (/) 


(4) 


where 


F^if)  =E{[n(f)+h{f)z{f)e~^^‘^^]  [A  (/)  +  ft  (/)  z  (/)  ’  } 


-2KiJTaJ  . 


=  F„(,f)+h(J)h^  (f)  .p(f) 


and  z  (/)  is  the  Fourier  transform  of  the  signal  waveform,  is  the  travel-time  to  the  refer¬ 
ence  sensor  of  the  array,  p(/)  is  the  power  spectral  density  of  the  signal  waveform  and  E 
denotes  statistical  expectation. 

Using  a  formula  of  Bartlett  ( 195 1 )  for  matrix  inversion,  we  can  write: 

F-J(f)  =  F;^{f)-F-^(f)h{f)h^  (f)F;^(f)p(f)[\+h^  *  (5) 

Substituting  (5)  into  (4)  yield.s: 

d(f)  =  - - - - 7 - ^  =  gU)  (6) 

(f)  [\-h^  {f)F;^(f)h(f)  -cj 

where  /  is  the  identity  matrix  and  c  =  p(j).[\+P  (/)F;'(/)ft(/)] 

Equation  (6)  shows  that  in  theory,  AOGF  will  work  even  if  the  adaptation  interval  con¬ 
tains  the  plane  wave  we  want  to  extract. 

Processing  of  synthetic  mixtures 

The  main  problem  arising  when  attempting  to  utilize  the  derived  theoretical  property  of 
AOGF  in  data  processing  practice,  is  to  estimate  F~'  (f)  with  sufficient  accuracy  using  a 
mixture  of  signal  and  noise  records.  Autoregressive  modelling  in  principle  gives  us  such 
an  opportunity.  This  will  be  shown  in  some  experiments  with  simulated  array  data. 

Using  the  geometry  of  the  NORESS  array,  we  have  synthesized  a  mixture  of  a  transient 
plane  wave  and  a  stationary  coherent  noise  field  with  some  added  white  noise.  The  slow¬ 
ness  vector  of  the  modelled  transient  is  identical  to  that  of  a  P-wave  originating  at  the 
Novaya  Zemlya  test  site.  With  relevance  to  the  problem  case  outlined  in  the  introduction, 
the  slowness  vector  of  the  stationary  coherent  noise  field  corresponds  to  that  of  a  P-wave 
from  an  event  in  the  Hindu  Kush  region.  Using  the  current  implementation  of  the  AOGF 
method,  a  successful  signal  extraction  is  shown  in  Fig.  7.5.  2.  The  conventional  beam 
shown  on  the  top  trace  did  not  reveal  the  signal,  whereas  the  signal  clearly  stands  out  on 
the  AOGF  filtered  trace  shown  in  the  middle.  For  comparison,  a  single  array  sensor  is 
shown  at  the  bottom.  This  experiment  shows  that  the  current  AOGF  implementation  is 
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capable  of  extracting  a  simple  signal  from  a  simple  noise  field,  even  if  adaptation  is  made 
to  the  data  containing  the  signal  we  want  to  extract. 

Experiments  with  NORESS  data  with  a  real  Novaya  Zemlya  signal  superimpo.sed  in  the 
coda  of  a  real  Hindu  Kush  event  did  not  give  the  desirable  result  with  the  current  imple¬ 
mentation  of  the  AOGF  method.  This  can  be  explained  by  the  specific  features  of  the  soft¬ 
ware,  which  has  been  designed  for  on-line  implementation,  as  shown  by  Kushnir  et  al. 
(1990).  In  this  program,  adaptation  is  made  after  steering  the  array  to  the  cho.sen  direction, 
resulting  in  synchronized  signal  waveforms  when  data  contain  the  signal  we  want  to 
extract.  Our  working  hypothesis  is  that  this  creates  a  problem  for  the  multichannel  autore¬ 
gressive  estimation  procedure,  in  particular  if  the  background  noise  has  rather  complex 
spatial  features  like  an  earthquake  coda. 

An  additional  experiment  indicates  that  if  we  introduce  the  appropriate  changes  to  the 
AOGF  implementation  (by  doing  all  computations  in  the  frequency  domain  without  steer¬ 
ing  of  the  array  sensors),  a  .successful  signal  extraction  can  be  achieved,  even  in  the  case 
when  the  adaptation  interval  contains  a  real  se'smic  signal  superimpo.sed  in  the  coda  of 
another  seismic  event. 

The  experiment,  simulating  the  frequency  domain  version  of  AOGF  by  use  of  the  current 
AOGF  software,  consi.sted  of  the  following  steps: 

1 .  A  synthetic  mixture  was  created  by  superimposing  down-.scaled  records  of  the 
Novaya  Zemlya  explosion  of  24  October  1990  in  the  coda  of  a  real  Hindu  Kush 
event  (event  1  of  Table  7.5.1).  The  SNR  between  the  Novaya  Zemlya  and  the  Hindu 
Kush  signals  was  approximately  1. 

2.  Adaptation  to  this  data  was  made  without  introducing  time  delays  (steering  with  infi¬ 
nite  apparent  velocity),  thus  avoiding  synchronized  waveforms. 

3.  A  new  .synthetic  mixture  was  created  by  superimposing  a  transient  plane  wave  with 
infinite  apparent  velocity  and  the  waveform  of  the  Novaya  Zemlya  explosion  in  the 
coda  of  the  Hindu  Kush  event.  In  this  case  the  SNR  was  only  0.25. 

4.  This  data  was  then  AOGF  filtered  with  the  parameters  obtained  during  the  previous 
adaptation. 

The  result  is  given  in  Fig.  7.5.3,  indicating  the  successful  extraction  of  the  signal. 

This  procedure  is  similar  to  what  will  be  done  by  the  new  frequency-domain  implementa¬ 
tion  of  AOGF,  thus  indicating  that  the  theoretical  property  outlined  in  the  preceding  sec¬ 
tion.  may  also  be  valid  in  practice. 

Adaptation  to  a  neighboring  event 

We  have  also  investigated  the  possibility  of  using  two  neighboring  events  in  the  AOGF 
proce.ssing,  using  one  event  for  adaptation,  and  subsequently  AOGF  filtering  the  other 
event  (which  may  contain  another  signal)  with  the  parameters  obtained  during  the  adapta¬ 
tion  to  the  first  event. 
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Two  Strong  Hindu  Hush  earthquakes  with  an  epicentral  difference  of  approximately  100 
km  were  chosen  for  this  experiment  (see  Table  7.5.1).  The  first  20  .seconds  of  the  two  P- 
wavetrains  recorded  at  NORESS  are  shown  in  Fig.  7.5.4,  and  the  spectra  of  the  first  10 
seconds  are  given  in  Fig.  7.5.5.  Note  that  the  spectra  show  very  similar  characteristics. 

The  results  of  this  test  are  summarized  in  the  spectra  of  Fig.  7.5.6.  The  upper  curve  shows 
the  spectrum  of  a  conventional  beam  steered  with  a  slowness  vector  corresponding  to  a  P- 
wave  from  the  Novaya  Zemlya  test  site,  using  NORESS  data  for  Hindu  Kush  event  no.  1 
(.see  Table  7.5.1 ).  The  spectrum  shown  by  the  solid  line  below,  corresponds  to  the  AOGF 
beam  steered  to  Novaya  Zemlya  of  Hindu  Kush  event  no.  1,  using  Hindu  Kush  event  no.  2 
for  adaptation;  and  the  spectrum  shown  by  the  dashed  line,  corresponds  to  the  AOGF 
beam  steered  to  Novaya  Zemlya  of  Hindu  Ku.sh  event  no.  1,  using  the  same  event  (no.  1) 
for  adaptation. 

Comparing  the  spectra  suggests  that  for  frequencies  below  1  Hz,  the  AOGF  method  can 
suppress  the  effect  of  interfering  events  much  better  than  conventional  beamforming, 
when  a  neighboring  event  has  been  used  for  adaptation.  For  higher  frequencies,  the 
improvement  is  also  significant.  The  spectrum  shown  by  the  dashed  line,  indicate  that 
more  improvement  can  be  obtained  if  the  same  event  is  used  for  adaptation. 

To  visualize  the  effect  of  using  a  neighboring  event  for  adaptation,  we  superimposed 
records  of  the  down-scaled  Novaya  Zemlya  explosion  in  the  coda  of  Hindu  Kush  event  no. 
1 .  As  in  Fig.  7.5.3,  the  SNR  was  0.25.  The  results  from  AOGF  filtering  this  mixture  signal 
with  the  adaptation  parameters  retrieved  from  Hindu  Kush  event  no.  2.,  are  shown  in  the 
middle  trace  of  Fig.  7.5.7.  As  expected,  e  to  signal  loss  during  the  plane- wave  beam¬ 
forming,  the  amplitude  of  the  Novaya  Zemlya  signal  is  reduced  with  comparison  to  the 
amplitude  of  Fig.  7.5.3.  Also  indicated  by  the  spectra  of  Fig.  7.5.6,  the  amplitude  of  the 
Hindu  Kush  event  is  not  as  effectively  suppressed  as  in  the  case  of  Fig.  7.5.3,  but  we  are 
still  able  to  identify  the  Novaya  Zemlya  arrival  on  the  ACX3F  beam. 

To  further  compare  the  relative  performance  between  the  conventional  and  the  AOGF 
beams,  we  bandpass  filtered  the  traces  of  Fig.  7.5.7  in  the  frequency  band  providing  the 
largest  SNR  (i.e.,  1. 5-3.5  Hz).  The  results  are  shown  in  Fig.  7.5.8.  The  only  difference 
from  Fig.  7.5.7,  is  that  instead  of  using  all  sen.sors  of  the  NORESS  array  when  forming  the 
conventional  beam,  we  used  a  sub-geometry  consisting  of  AOZ,  the  C-ring  and  the  D-ring 
(17  sensors).  In  the  frequency  band  1. 5-3.5  Hz,  this  sub-geometry  will  provide  the  best 
SNR  improvement,  at  lea.st  when  comparing  to  background  noise  conditions  (Kvaema, 
1989).  The  result  shows  that  when  using  a  neighboring  event  for  adaptation,  the  AOGF 
beam  provides  better  SNR  than  the  optimum  conventional  beam  (optimum  =  best  sub¬ 
geometry  and  best  filter  band). 

Discussion 

The  results  presented  above  suggest  that  application  of  the  A(X3F  method  has  a  potential 
for  significantly  improving  the  signal-to-noise  ratios  of  seismic  events  occurring  in  the 
coda  of  large  earthquakes.  This  is  important,  e.g.,  for  improving  the  threshold  monitoring 
capability  during  periods  of  increa.sed  magnitude  thresholds.  The  experiments  with  the 
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Hindu  Kush  events  recorded  at  NORESS  indicate  that  with  the  appropriate  changes  to  the 
AOGF  software  (by  implementing  the  frequency  domain  version),  :i  significant  threshold 
magnitude  reduction  may  be  achieved  for  each  phase  considered. 

Applying  the  software  currently  available,  using  adaptation  to  a  neighboring  event,  we 
have  the  capability  of  some  reduction  of  the  magnitude  thresholds. 

The  continuation  of  this  research  will  consist  of  caiefully  testing  the  frequency-domain 
version  of  ACXjF,  and  special  emphasis  will  be  devoted  to  ensuring  the  undistorting  fea¬ 
ture  of  AOGF  on  the  signal  we  want  to  extract  from  the  background  coda.  In  addition, 
practical  problems  related  to  the  use  of  AOGF  in  the  threshold  monitoring  algorithm  will 
be  investigated. 


Tormod  Kvxrna 

Alex  Kushnir,  MITPAN  Institute,  Moscow,  USSR 
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Event  no. 

Year 

Day 

Mon 

Ortime 

Lat 

Lon 

Depth 

•"b 

1 

1990 

2S 

Oct 

04.53.59.9 

35.121N 

70.486E 

114 

6.0 

2 

1990 

15 

May 

14,25.20.6 

36.043N 

70.428E 

113 

5.9 

Table  7,5.1.  PDE  locations  of  the  Hindu  Kush  events. 
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Fig.  7.5.1.  Threshold  monitoring  of  the  Novaya  Zcmiya  test  site  for  day  298(25  October  1990). 

The  top  three  traces  represent  thresholds  (upper  90  per  cent  magnitude  limits)  obtained  from 
each  of  the  three  arrays  (ARCESS,  FINESA,  NORESS),  whereas  the  bottom  trace  shows  the 
combined  network  thresholds.  Tlic  FINESA  array  had  several  short  outages  this  day,  but  this 
caused  no  particular  problems  in  tenns  of  network  threshold  capacity. 

Notes: 

3.  An  earthquake  (m(,  =  4.5)  near  Jan  Mayen.  The  corresponding  network  threshold  peak  for 
Novaya  Zcmiya  is  mj,  =  2.8. 

4.  A  tclcscismic  earthquake  (m^  =  6.0)  at  Hindu  Kush.  The  relatively  strong  P-wavc  train 
caused  a  peak  threshold  of  mj,  =  3.8  for  monitoring  Novaya  Zcmiya. 

5.  A  tclcscismic  earthquake  (nif,  =  5.9)  at  Mindanao.  Philippine  Islands.  Corresponding  thresh¬ 
old  is  mf,  =  3.0.  i 

6. -7.  A  sequence  of  seismic  events  (presumably  underwater  explosions)  near  Mumiansk,  Kola  I 

Peninsula.  The  network  threshold  for  monitoring  Novaya  Zcmiya  is  about  rri),  =  2.5  to  2.8  at  [ 

the  times  of  these  events.  ; 
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Fig.  7.5.2.  Synthetic  processing  example. 

A  synthetic  signal  has  been  created  by  adding  three  components  as  follows: 

1 )  A  stationary,  coherent  noise  field  (Hindu  Kush  direction,  P  velocity). 

2)  A  stationary,  incoherent  "white  noise”  field. 

3)  A  transient,  coherent  signal  (Novaya  Zemlya  direction,  P  velocity). 

The  ton  trace  shows  the  conventional  beam  steered  to  Novaya  Zemlya. 
Thejniddle  trace  shows  the  ACXrF  beam  steered  to  Novaya  Zemlya. 

The  bottom  trace  shows  the  single  sensor  AOZ. 

Note  that  the  signal  is  clearly  visible  on  the  AOGF  trace,  but  not  on  the  other  two. 
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Fig.  7S3.  Synthetic  processing  example,  illustrating  that  with  the  appropriate  changes  to  the 
AOGF  program,  we  may  have  the  possibility  of  extracting  a  real  seismic  signal  from  an 
earthquake  coda,  even  if  the  adaptation  interval  containing  the  signal  we  want  to  extract 

The  experiment,  simulating  the  frequency  domain  implementation  of  ACXjF,  consisted  of  the  fol¬ 
lowing  steps: 

1 )  A  synthetic  mixture  was  created  by  superimposing  down-scaled  records  of  the  Novaya  Zemlya 

explosion  of  24  October  1 990  in  the  coda  of  a  real  Hindu  Kush  event  (event  1  of  Table  7.5.1). 
The  SNR  between  the  Novaya  Zemlya  and  the  Hindu  Kush  signals  was  approximately  1 . 

2)  Adaptation  to  this  data  was  made  without  introducing  time  delays  (steering  with  infinite  appar¬ 

ent  velocity),  thus  avoiding  synchronized  waveforms. 

3)  A  new  synthetic  mixture  was  created  by  superimposing  a  transient  plane  wave  with  infinite 

apparent  velocity  and  the  waveform  of  the  Novaya  Zemlya  explosion  in  the  coda  of  the 
Hindu  Ku.sh  event.  In  this  case  the  SNR  was  only  0.25. 

4)  This  data  was  then  A(XJF  filtered  with  the  parameters  obtained  during  the  previous  adaptation. 

The  top  trace  shows  the  conventional  beam  .steered  with  infinite  apparent  velocity. 

The  middle  trace  shows  the  AOGF  beam  steered  with  infinite  apparent  velocity. 

The  bottom  trace  shows  the  single  sensor  AOZ. 

Note  that  the  signal  is  clearly  visible  on  the  AOGF  trace,  but  not  on  the  other  two. 
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Fig.  7^.4.  NORESS  recordings  (instrument  AOZ)  of  the  two  Hindu  Kush  events  described  in 
Table  7.5.1.  Event  no.  1  is  shown  in  the  upper  trace,  and  event  no.  2  in  the  lower. 
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NORESS  recordings  of  the  two  Hindu  Kash  events  described  in 

Table  7^5. ! .  SqM  IiDC:  Event  no.  I  Dastied  line:  Event  no.  2.  Both  spectra  were  estimated 
using  1 0  seconds  of  data. 
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Fig.  7.5.6.  After  processing  Hindu  Kush  event  no.  1  with  different  techniques,  we  computed  P- 
wave  amplitude  spectra  of  the  first  10  seconds  of  the  signal.  The  unper  spectrum  results 
from  conventional  beamforming  using  the  Novaya  Zemlya  P-wave  steering  delays.  The 
loyCf  solid  line  ^pwprutn  rosults  from  AOGF  filtering  with  Novaya  Zemlya  P-wave  steer¬ 
ing  delays  using  Hirtdu  Kush  event  no.  2  for  adaptation.  The  dashed  line  spectrum  rrmihg 

from  ACXjF  filtering  with  Novaya  Zemlya  P-wave  steering  delays  using  the  same  event  (no. 
1 )  for  adaptation. 
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Fig.  IS.l.  Results  after  processing  the  mixture  of  Hindu  Kush  event  no.  1  and  the  down-scaled 
Novaya  Zemlya  signal.  Single  channel  SNR  ~  0.25.  The  lower  trace  shows  a  single 
NORESS  sensor,  the  top  trace  is  a  convcntiOTial  beam  steered  to  Novaya  Zemlya,  whereas 
the  middle  trace  is  the  ACXjiF  output  using  Hindu  Kush  event  no.  2  for  adaptation. 
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60,000  Sec 


Fig.  7.5.8.  Results  from  bandpass  filtering  the  data  of  Fig.  7.5.7  in  the  passband  1 .5-3.5  Hz.  The 
only  difference  from  Fig.  7.5.7,  is  that  instead  of  using  all  sensors  of  the  NORESS  array  when 
forming  the  conventional  beam,  we  used  a  sub-geometry  consisting  of  AOZ,  the  C-ring  and 
the  D-ring  ( 1 7  sensors).  Note  Uiat  even  under  these  optimum  conditions  for  cmventional 
beamforming,  the  AOGF  beam  provide  the  best  SNR. 
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7.6  Regional  detection  performance  during  GSETT-2:  Initial  results 
for  the  Fennoscandian  array  network 


Introduction 

During  the  period  22  April  to  9  June  1991,  the  Conference  on  Disarmament’s  Group  of 
Scientific  Experts  carried  out  the  main  phase  of  its  Second  Technical  Test  (GSETT-2).  A 
total  of  34  countries  participated  in  this  test,  providing  seismic  data  for  42  consecutive 
data  days  from  60  stations  di.stributed  around  the  globe  (Fig.  7.6.1).  Data  were  recorded 
and  proces.sed  at  National  Data  Centers,  and  results  as  well  as  waveform  segments  were 
transmitted  to  four  experimental  International  Data  Centers  (EIDCs)  for  further  analysis. 
Results  of  these  analyses  were  summarized  in  event  bulletins,  which  were  transmitted 
back  to  participants  from  the  EIDCs. 

An  important  aspect  of  the  performance  evaluation  of  GSETT-2  is  the  completeness  and 
quality  of  the  final  event  bulletin  (FEB).  This  seismological  output  is  clo.sely  linked  to  the 
actual  spatial  distribution  of  .seismic  stations.  For  GSETT-2,  a  very  heterogeneous  global 
coverage  yielded  large  regional  variations  in  detection  threshold.  About  one  half  of  the 
participating  stations  were  situated  in  and  around  Europe,  consequently  a  lai^e  number  of 
small  events  were  detected,  tly  quarry  blasts  and  rock  bursts  of  magnitude  1  to  4.  On 
the  other  hand,  in  many  areas  of  the  globe  where  the  .station  distribution  was  very  sparse, 
only  latter  earthquake.s  ere  detected. 

In  this  paper  a  preliminary  assessment  is  made  of  the  detection  capability  during  GSETT-2 
for  Fennoscandia  and  NW  Russia.  This  is  the  region  that  had  maybe  the  best  instrumental 
coverage  during  the  experiment.  In  particular  the  regional  arrays  deployed  in  this  area 
made  significant  contributions. 

Our  results  for  this  region  represent  in  a  sense  the  “best"  regional  performance  during 
GSETT-2.  It  is  in  no  way  representative  for  the  performance  on  a  global  or  more  extended 
regional  scale.  However,  it  does  serve  to  illu.strate  the  potential  capabilities  of  a  monitor¬ 
ing  network,  assuming  that  an  adequate  density  and  number  of  high  quality,  sensitive 
array  stations  are  deployed. 

Method 

The  method  u.sed  for  detectability  estimation  has  been  described  by  Ringdal  (1975),  and  is 
briefly  summarized  as  follows: 

1 .  A  reference  system,  independent  of  the  system  to  be  evaluated,  is  used.  Event  lists 
and  magnitudes  from  this  reference  system  are  compiled. 

2.  For  each  reference  event,  a  comparison  is  made  to  see  if  the  system  to  be  evaluated 
has  detected  the  event, 

3.  Ba.sed  on  the  number  of  detections/no  detections  at  each  magnitude,  a  maximum 
likelihood  approach  is  made  to  estimate  a  “detection  curve"  of  the  form 
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^  e  dx  (1) 

Jlna 

Here  G  (/w;^,  o)  denotes  the  incremental  probability  of  detection,  given  event 
magnitude  m.  The  detection  curve  is  completely  characterized  by  the  parameters  p. 
and  a.  The  50  and  90  per  cent  incremental  detection  thresholds  and  p^) 
become; 


G(m;p,  a)  =  | 


1^50  = 

(2) 

=  p+ 1.29-0 

(3) 

It  should  be  noted  that  while  the  method  assumes  that  the  reference  network  pro¬ 
vides  independent  event  estimates,  it  is  not  necessary  to  have  a  complete  event  cat¬ 
alogue  in  any  given  magnitude  range.  Thus  the  reference  events  actually  selected 
are  assumed  to  be  randomly  sampled  from  the  total  number  of  events  available, 
much  in  the  same  way  as  opinion  survey  polls  attempt  to  address  randomly 
selected  subsets  of  the  population.  The  resulting  detectability  estimates  will  be  rep¬ 
resentative  for  the  region  considered  only  to  the  extent  that  the  reference  event  set 
is  representative. 

Reference  network 

The  reference  data  base  for  this  study  has  been  the  catalogue  of  seismic  events  in  northern 
Europe  regularly  compiled  by  the  Seismological  Institute,  University  of  Helsinki. 

The  stations  used  in  compiling  this  catalogue  are  in  almost  all  cases  comprised  of  the 
Finnish  seismic  network  single  .stations.  For  all  practical  purposes,  the  compilation  is 
independent  of  the  regional  arrays  in  Fennoscandia  (NORESS,  ARCESS,  FINESA).  The 
magnitudes  quoted  in  the  bulletin  are  likewise  derived  independently  of  the  regional 
arrays,  and  comprise  either  duration  magnitudes  (in  most  cases)  or  local  magnitudes. 
These  magnitudes  are  fairly  consistent  with  magnitudes  calculated  by  the  Intelligent  Mon¬ 
itoring  System,  while  their  relationship  to  teleseismic  m^,  estimates  is  at  present  not  well 
established. 

For  the  month  of  May  1991,  upon  which  this  analysis  is  based,  the  reference  catalogue 
contained  321  seismic  events  in  the  region  bounded  by  58°-70°N,  20°-40®E,  of  which  108 
had  an  assigned  magnitude  in  the  range  1. 7-2.9. 

Results 

The  initial  results  from  the  detectability  study  are  presented  in  Figs.  7.6.2-7.6.4.  Each  fig¬ 
ure  is  based  upon  analyst  comparison  of  the  reference  events  with  bulletin  reports  accord¬ 
ing  to  different  reporting  criteria; 


I2H 


NORSAR  Sci.  Rep.  1-91/92 


November  1991 


a)  One-arrav  detection: 

In  Fig.  7.6.2,  an  event  is  considered  detected  if  it  was  reported  with  2  phases  (P  and 
S;  or  P  and  Lg)  by  at  least  one  of  the  three  regional  arrays  (NORESS,  ARCESS, 
FINESA).  In  terms  of  GSETT-2  final  event  bulletins,  this  means  that  the  event 
would  either  be  located  as  a  multi-station  event,  or  listed  as  an  NDC-reported 
event.  We  note  that  the  50  %  threshold  is  close  to  1.7,  and  the  90  %  threshold  is  2.3 
in  this  case. 

b)  Two-arrav  detection: 

In  Fig.  7.6.3,  the  detection  requirement  is  2  pha.ses  (P  and  S;  or  P  and  Lg)  from  one 
array  and  at  least  one  confirming  phase  from  another  array.  This  added  requirement 
has  the  effect  of  increasing  the  50  %  threshold  to  2.1 .  and  the  90  %  threshold  to  2.4. 

c)  GSETT-2  reported  events: 

In  Fig.  7.6.4,  the  GSETT-2  reported  events,  located  by  at  least  one  IDC,  are  shown. 
(We  have  not  counted  as  detected  events  those  events  whose  definition  depended 
upon  reportings  from  the  Finnish  network  stations  KAF  and  VAF,  since  these  two 
stations  were  part  of  the  reference  network.)  The  resulting  thresholds  are  similar  to 
those  displayed  in  Fig.  7.6.3  (the  two-array  case). 

Conclusions 

The  regional  evaluation  of  detection  results  from  GSETT-2  presented  here  shows  that  in  a 
region  with  dense  coverage  of  high-quality  array  stations  as  in  Fennoscandia,  it  is  possible 
to  detect  seismic  events  at  very  low  magnitudes. 

The  90  per  cent  threshold  of  around  magnitude  2.5  found  in  this  study  must  of  course  be 
considered  with  the  appropriate  caution;  thus  it  refers  to  a  regional  magnitude  scale  that 
currently  is  not  well  calibrated  in  terms  of  global  magnitude.  Also,  in  other  geological 
environments,  the  wave  propagation  and  array  noi.se  suppression  characteristics  may  be 
different.  Therefore,  it  is  not  known  to  which  extent  such  results  would  be  possible  to 
duplicate  in  other  parts  of  the  world. 

Further  work  should  be  undertaken  to  evaluate  the  GSETT-2  detection  performance  in 
other  regions,  in  order  to  gain  experience  for  future  studies.  It  will  also  be  necessary  to 
study  event  location  performance,  and  evaluate  how  the  detection  capability  varies  as  a 
function  of  different  requirements  with  regard  to  event  location  accuracy. 

F.  Ringdal 
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Fig.  7.6.1.  f  IDCs.  NDCs  and  stations  participating  in  GSETT-2,  April  June  1991 
(After  GSEAJS/68). 
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Fig.  7.6.2.  Maximum  likelihood  detecf  'Hility  estimation  for  Fennoscandia-NW  Russia 

using  the  Univ.  of  Helsinki  bulletin  a'^  ference.  The  upper  half  shows  the  reference 

event  set  and  the  number  of  events  actually  detected  for  each  magnitude.  The  lower  half 

shows  the  maximum  likelihood  detectability  curve  and  its  confidence  limits.  The  actual  | 

percentage  of  detected  events  at  each  magnitude  is  also  shown.  This  figure  is  based  upon  a 

one-array  detection  requirement. 
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Fig.  7.6J.  Maximum  likelihood  detectability  estimation  for  Fennoscandia-NW  Russia 
using  the  Univ.  of  Helsinki  bulletin  as  a  reference.  The  upper  half  shows  the  reference 
event  set  and  the  number  of  events  actually  detected  for  each  magnitude.  The  lower  half 
shows  the  maximum  likelihood  detectability  curve  and  its  confidence  limits.  The  actual 
percentage  of  detected  events  at  each  magnitude  is  also  shown.  This  figure  is  based  upon  a 
two-array  detection  requirement. 
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F^.  7.6.4.  Maximum  likelihood  detectability  estimation  for  Fennoscandia-NW  Russia 
using  the  Univ.  of  Helsinki  bulletin  as  a  reference.  The  upper  half  shows  the  reference 
event  set  and  the  number  of  events  actually  detected  for  each  magnitude.  The  lower  half 
shows  the  maximum  likelihood  detectability  curve  and  its  confidence  limits.  The  actual 
percentage  of  detected  events  at  each  magnitude  is  also  shown.  This  figure  is  based  upon 
FEB  reported  events  as  discussed  in  the  text. 
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7.7  Contemporary  seismicity  of  the  N  W  part  of  the  USSR 

Introduction 

The  regional  seismic  network  of  the  Kola  Science  Centre  began  to  be  established  in  1982. 
Today,  7  seismic  stations  are  used  to  monitor  the  seismic  activity  of  the  Barents  Sea  and 
the  NW  part  of  the  USSR.  The  station  list  is  shown  in  Table  7.7.1.  All  stations  are 
equipped  with  short-period  seismometers  with  photo  sensitive  analog  recording,  and  paper 
speed  60  mm/min.  A  three-component  digitally  recording  station  has  recently  (June  1991) 
been  installed  in  Apatity.  For  better  location  of  seismic  events  in  the  Barents  Sea  platform, 
data  from  the  seismic  station  KHE  in  Franz  Josef  Land  (80.37®N,  58.03°E)  are  used,  even 
though  this  station  is  not  part  of  the  network. 

Strong  earthquakes  are  known  from  as  far  back  as  the  17th  century  for  the  eastern  Baltic 
Shield  (T.  Ahjos  and  H.  Korhonen,  1984),  and  before  the  seismic  networks  in  this  region 
were  installed,  it  was  suggested  that  only  large  earthquakes  could  take  place  here,  because 
the  crust  of  the  rigid  ancient  shield  was  believed  not  to  be  able  to  produce  many  small 
earthquakes  (e.g.,  Panasenko,  1969).  Now  this  opinion  has  been  revised. 

The  contemporary  seismicity  in  northern  Europe  is  shown  in  Fig.  7.7.1.  Most  of  the  earth¬ 
quakes  shown  are  for  the  past  30-year  period,  but  some  historical  epicenters  are  also 
included,  dating  back  as  far  as  the  year  1542.  The  earthquake  activity  is  characterized  by 
small  and  moderate  size  earthquakes.  The  largest  earthquake  in  the  Baltic  Shield  in  recent 
years  was  the  mL  =  5.2  event  of  20  May  1%7,  located  at  66.6°N,  33.7®E  (Meyer  and 
Ahjos,  1985).  It  was  felt  with  MSK  intensity  VII  in  Karelia,  and  generated  a  small  tsunami 
in  the  White  Sea. 

The  main  seismic  zones  of  the  Kola  Peninsula  are  seen  more  clearly  when  considering 
only  the  time  interval  1986  to  1991.  They  comprise  the  Kola-Finnmark  zone  and  the  Khib- 
iny  Massif  (Fig.  7.7.2). 

The  Kola-Finnmark  zone 

Information  about  the  seismicity  of  the  Kola-Finnmark  zone  was  known  already  from 
1772  when  an  earthquake  with  maximum  intensity  VI  took  place  near  Murmansk  (68.7°, 
33.3°E).  The  laig:est  earthquake  known  in  the  Kola-Finnmaik  zone  (Kola  Peninsula)  dur¬ 
ing  this  century  occurred  on  10  April  1981  (68.8°N,  36.0°E).  The  magnitude  estimated  for 
this  event  ranges  from  4.5  Ml  (Upp)  to  4.7  mb  (ISC).  The  1981  earthquake  had  two 
shocks  about  10  seconds  apart  and  was  felt  with  intensity  V  at  Murmansk. 

During  1986-1990  some  tens  of  events  were  registered  from  the  Kola-Finnmaiic  zone.  One 
of  the  larger  of  these  was  the  earthquake  that  occurred  on  16  June  1990  (69.14°N, 
35.15°E)  with  Ml  =  4.0,  which  was  felt  with  intensity  IV  at  Murmansk. 

The  seismicity  of  the  Murmansk-Finnmark  zone  is  laigely  confined  to  intersections  of 
faults  separating  the  Murmansk  Block  from  the  South  Barents  Sea  platform,  i.e.,  the  so- 
called  Karpinsky  lineament,  and  to  NS-trending  faults  which  are  clearly  discerned  from 
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Landsat  imagery  (Fig.  1.1. 'i).  Fig.  7.7.4  presents  a  fragment  of  the  neotectonic  framework 
for  these  earthquakes. 

The  Khibiny  Massif 

In  the  Ia.st  several  years  there  has  been  a  great  increase  of  induced  seismicity  in  the  Khib¬ 
iny  Massif.  This  began  to  manifest  itself  after  lat^e  volumes  of  rocks  had  been  extracted 
from  a  high-.stress  environment.  Fig.  7.7.5  demonstrates  the  position  of  the  five  main 
mines  of  the  Khibiny  Massif  (a)  and  allows  a  comparison  of  the  quantity  of  rock  removal 
with  the  number  of  Khibiny  earthquakes  per  year  (b). 

Exploration  of  the  Khibiny  apatite  deposits  was  started  in  1929.  At  the  present  time  the 
production  is  more  than  50  mill,  tons  per  year.  The  first  earthquake  on  record  in  the  Khib¬ 
iny  area  occurred  on  23  September  1948  and  was  felt  with  intensity  IV  in  populated  areas. 

The  largest  earthquake  in  Khibiny  took  place  on  16  April  1989  (Ml  =  4.2).  This  event  was 
felt  with  intensity  VIII  in  the  upper  levels  of  the  Kirovsk  mine  and  caused  considerable 
destruction.  The  maximum  measured  displacement  was  15-20  cm  and  it  occurred  along  a 
fault  striking  at  125-135°  and  dipping  at  30-35°  NE.  This  displacement  was  traced  along 
the  surface  for  1200  m  and  observed  to  a  depth  of  at  least  220  m. 

The  earthquake  was  felt  with  intensity  V-VI  in  the  Kirovsk  area,  and  it  was  followed  by 
several  hundred  aftershocks  during  the  next  two  months.  The  earthquake  on  16  April  1989 
occurred  simultaneously  with  a  240  ton  explosion  in  one  of  the  Kirovsk  mines.  It  thus 
appears  that  the  explosion  triggered  the  earthquake.  Supporting  this  assertion  is  the  obser¬ 
vation  of  analogous  events  at  the  Kirov.sk  mine  on  29  August  1982,  when  an  earthquake 
with  maximum  intensity  VI  took  place  at  the  same  time  as  a  large  explosion.  A  similar 
pair  of  events  occurred  on  28  March  1991. 

The  location  of  mines  and  earthquakes  relative  to  each  other  is  shown  in  Fig.  7.7.6. 

NE  Archangelsk  and  Novaya  Zemlya 

The  .seismic  station  Amderma  (AMD)  began  operation  in  1983  (see  Fig.  7.7.7).  The  sta¬ 
tion  was  established  to  study  seismic  activity  of  the  Barents  Sea  Platform.  The  first  inves¬ 
tigations  identified  two  zones  of  weak  seismic  activity  in  the  NE  Archangelsk  region 
(Assinovskaya,  1989).  The  shaded  area  in  Fig.  7.7.7  shows  the  .seismic  zone  to  the  south¬ 
east  of  Amderma.  No  definite  information  exists  about  the  nature  of  the  events  recorded  in 
this  region,  but  neotectonic  conditions  are  consi.stent  with  this  zone  being  .seismically 
active. 

A  .seismically  active  zone  has  also  been  found  on  Novaya  Zemlya.  Of  particular  interest  is 
the  earthquake  on  1  August  1 986  with  m^,  =  4.6  near  the  south  coast  of  the  Matochkin  Shar 
strait.  The  location  of  this  event  is  shown  in  Fig.  7.7.8  together  with  the  system  of  linea¬ 
ments  and  faults.  There  were  at  least  four  more  earthquakes  in  this  area,  information  about 
which  has  been  included  in  Table  7,7.2. 
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In  conclusion,  we  find  that  the  seismicity  of  the  NW  part  of  the  USSR  is  characterized  by 
clustering  of  the  most  intense  earthquakes  at  sites  of  intersections  of  NS  and  EW  trending 
fault  zones. 


E.O.  Kremenctskaya,  Kola  Science  Centre,  USSR 
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Name 

Coordinates 

1 

APA 

67.568N 

33.388E 

AMD 

69.742N 

61.655E 

- 

BRB 

78.073N 

14.240E 

PLZ 

67.400N 

32.533E 

PLQ 

66.4  ION 

32.750E 

PLG 

62.320N 

36.930E 

KEM 

64.956N 

34.635E 

Table  7.7.1.  List  of  seismic  stations  operated  by  the  Kola  Science  Centre. 


Date  Time  Coordinate  m|, 

N  E 


20.07.83 

015941.7 

73.2 

56.3 

- 

31.01.86 

110813.0 

73.2 

56.3 

2.1 

21.04.86 

093753.9 

73.2 

56.3 

3.2 

01.08.86 

135636.8 

73.18 

56.32 

4.6 

29.12.86 

034407.4 

- 

- 

3.0 

Table  7.7.2.  Recent  earthquakes  in  Novaya  Zemlya. 
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Fig.  7.7.1 .  Epicentral  map  of  earthquakes  in  the  European  Arctic  region.  Most  of  the  earth¬ 
quakes  shown  are  from  the  period  1960-1990,  but  some  earlier  events  have  been  included. 
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Fig.  7.7.2.  Epicentral  map  of  earthquakes  of  the  eastern  part  of  the  Baltic  Shield  1986- 
1991. 
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Fig.  7.7  Epicentral  map  of  earthquakes  of  the  Barents  Sea  with  tectonic  features: 

1)  faults  from  geophysical  data;  2)  faults  from  geological  data;  3)  thrusts;  4)  dislocations. 
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Fig.  7.7.4.  Neotectonic  framework  for  earthquakes  of  the  Murmansk  coast:  1)  faults  from 
geophysical  data;  2)  dislocations;  3)  deep  structures.  The  epicenter  of  the  earthquake  of  10 
April  1981  discussed  in  the  text  is  marked  a.s  a  circle  near  the  center  of  the  map. 
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33.00  33.25  33-50  33.75  34.00  34.25  34.50 


Fig.  7.7.6.  Relative  position  of  earthquakes  and  mines  in  the  Khibiny  Massif.  The  location 
of  the  Apatity  seismic  station  (APA)  is  also  shown.  The  mines  (I-V)  are  shown  in  more 
detail  in  Fig.  7.7.5. 
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Fig.  7.7.8.  The  position  of  the  earthquake  of  01.08.86  with  some  neotectonic  features: 
1 )  faults;  2)  boundaries  of  blocks;  3)  epicenter. 


